Multiplexed digital assay with data exclusion for calculation of target levels

ABSTRACT

Digital assay system, including methods and apparatus, for calculating a level of one or more targets. In an exemplary method, data for amplification of a plurality of targets including a first target may be collected from partitions. A level of the first target may be calculated from only a subset of the data that excludes partitions according to target content.

CROSS-REFERENCES TO PRIORITY APPLICATIONS

This application is based upon and claims the benefit under 35 U.S.C.§119(e) of U.S. Provisional Patent Application Ser. No. 61/759,772,filed Feb. 1, 2013; U.S. Provisional Patent Application Ser. No.61/759,930, filed Feb. 1, 2013; and U.S. Provisional Patent ApplicationSer. No. 61/759,931, filed Feb. 1, 2013. Each of these priorityapplications is incorporated herein by reference in its entirety for allpurposes.

CROSS-REFERENCES TO OTHER MATERIALS

This application incorporates by reference in their entireties for allpurposes the following materials: U.S. Pat. No. 7,041,481, issued May 9,2006; U.S. Patent Application Publication No. 2010/0173394 A1, publishedJul. 8, 2010; U.S. Patent Application Publication No. 2011/0217712 A1,published Sep. 8, 2011; U.S. Patent Application Publication No.2012/0152369 A1, published Jun. 21, 2012; U.S. Patent ApplicationPublication No. 2013/0040841 A1, published Feb. 14, 2013; U.S. patentapplication Ser. No. 14/099,750, filed Dec. 6, 2013; and Joseph R.Lakowicz, PRINCIPLES OF FLUORESCENCE SPECTROSCOPY (2^(nd) Ed. 1999).

INTRODUCTION

Digital assays generally rely on the ability to detect the presence oractivity of individual copies of an analyte in a sample. In an exemplarydigital assay, a sample is separated into a set of partitions, generallyof equal volume, with each containing, on average, less than about onecopy of the analyte. If the copies of the analyte are distributedrandomly among the partitions, some partitions should contain no copies,others only one copy, and, if the number of partitions is large enough,still others should contain two copies, three copies, and even highernumbers of copies. The probability of finding exactly 0, 1, 2, 3, ormore copies in a partition, based on a given average concentration ofanalyte in the partitions, is described by a Poisson distribution.Conversely, the concentration of analyte in the partitions (and thus inthe sample) may be estimated from the probability of finding a givennumber of copies in a partition.

Estimates of the probability of finding no copies and of finding one ormore copies may be measured in the digital assay. Each partition can betested to determine whether the partition is a positive partition thatcontains at least one copy of the analyte, or is a negative partitionthat contains no copies of the analyte. The probability of finding nocopies in a partition can be approximated by the fraction of partitionstested that are negative (the “negative fraction”), and the probabilityof finding at least one copy by the fraction of partitions tested thatare positive (the “positive fraction”). The positive fraction or thenegative fraction then may be utilized to determine the concentration ofthe analyte in the partitions, such as with Poisson statistics.

Digital assays frequently rely on amplification of a nucleic acid targetin partitions to enable detection of a single copy of an analyte.Amplification may be conducted via the polymerase chain reaction (PCR),to achieve a digital PCR assay. The target amplified may be the analyteitself or a surrogate for the analyte generated before or afterformation of the partitions. Amplification of the target can be detectedoptically from a fluorescent probe included in the reaction. Inparticular, the probe can include a fluorophore that provides afluorescence signal indicating whether or not the target has beenamplified.

A digital PCR assay can be multiplexed to permit detection of two ormore different targets within each partition. Amplification of thetargets can be distinguished by utilizing target-specific probes.However, in multiplexed assays with a higher concentration of targets,where two or more targets can be present together in the same partition,distinct partition populations within the assay often produceoverlapping signal amplitudes. Data interpretation for these assays canbe difficult, with calculated target concentrations being inaccurate.

New approaches are needed for improving multiplexed digital assays, suchas assays in which the data are not well resolved.

SUMMARY

The present disclosure provides a multiplexed digital assay system,including methods and apparatus, for calculating a level of one or moretargets. In an exemplary method, data for amplification of a pluralityof targets including a first target may be collected from partitions. Alevel of the first target may be calculated from only a subset of thedata that excludes partitions according to target content.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a flowchart of an exemplary method of performing a multiplexeddigital assay on targets in partitions containing a generic reporter,optionally with at least one of the targets being an obscuring targetthat at least partially masks the presence of another target inpartitions containing both the obscuring target and the masked target,in accordance with aspects of the present disclosure.

FIG. 2 is a schematic view of an exemplary system for performing themultiplexed digital assay of FIG. 1, in accordance with aspects of thepresent disclosure.

FIG. 3 is a schematic graph of multiplexed assay data that may becollected in a single detection channel from partitions, such asdroplets, collectively containing an obscuring target (H) associatedwith a higher amplitude signal and a masked target (L) associated with alower amplitude signal, in accordance with aspects of the presentdisclosure.

FIG. 3A is a schematic graph of multiplexed assay data that may becollected as in FIG. 3, but with the population of double-positivepartitions (HL, HLL, etc.) exhibiting a small, detectable increase insignal amplitude over the population of partitions containing onlytarget H, with the increase in signal amplitude being insufficient toresolve the two higher amplitude populations from each other, such thattarget L is still at least partially masked by the presence of target H,in accordance with aspects of the present disclosure.

FIG. 4 is another schematic graph of multiplexed assay data that may becollected generally as in FIG. 3, but with the obscuring target (L)having a lower amplitude signal than the masked target (H), inaccordance with aspects of the present disclosure.

FIG. 5 is a schematic graph of multiplexed assay data that may becollected in a single detection channel from partitions, such asdroplets, collectively containing K targets that mask one anotherserially according to a masking hierarchy, in accordance with aspects ofthe present disclosure.

FIG. 6 is a schematic view of a partition from an exemplary multiplexeddigital assay capable of generating the data of FIG. 3, with thepartition containing a copy of a masking target (H) and a copy of amasked target (L) of different length, each amplifiable at a differentefficiency with the same pair of forward (F) and reverse (R) primers anddetectable via the same probe (P), in accordance with aspects of thepresent disclosure.

FIG. 7 is a schematic view a partition from an exemplary multiplexeddigital assay capable of generating the data of FIG. 4, with thepartition containing a copy of a masking target (L) and a masked target(H) of different length, each amplifiable at a different efficiency withthe same pair of forward (F) and reverse (R) primers and detectable viathe same generic reporter for amplicons (e.g., an intercalating dye), inaccordance with aspects of the present disclosure.

FIG. 8 is a schematic view of a partition from another exemplarymultiplexed digital assay capable of generating the data of FIG. 3, withthe partition containing a copy of a masking target (H) and a copy of amasked target (L) of the same length, each amplifiable at a differentefficiency with the same pair of forward (F) and reverse (R) primers anddetectable via the same probe (P), in accordance with aspects of thepresent disclosure.

FIG. 9 is a schematic view of a partition from yet another exemplarymultiplexed digital assay capable of generating the data of FIG. 3, withthe partition containing a copy of a masking target (H) and a copy of amasked target (L) of the same length, each amplifiable at the sameefficiency with the same pair of forward (F) and reverse (R) primers anddetectable at different efficiencies via the same probe (P), inaccordance with aspects of the present disclosure.

FIG. 10 is a schematic view of a partition from still yet anotherexemplary multiplexed digital assay capable of generating the data ofFIG. 3, with the partition containing a copy of an obscuring target anda copy of a masked target, each amplifiable at a different efficiencywith a different pair of forward and reverse primers (F_(H) and R_(H) orF_(L) and R_(L)) and detectable via different probes (P_(H) and P_(L)),in accordance with aspects of the present disclosure.

FIG. 11 is a schematic view of a partition from an exemplary multiplexeddigital assay performed with a single detection channel, with thepartition containing a copy of target A and target B, each amplifiablewith a different pair of forward and reverse primers (F_(A) and R_(A) orF_(B) and R_(B)) and detectable via a generic reporter (e.g., anintercalating dye), in accordance with aspects of the presentdisclosure.

FIG. 11A is a schematic graph of fluorescence amplitude data that may becollected from five sets of partitions (lanes 1-5) in a single detectionchannel for the multiplexed digital assay of FIG. 11, with each setcontaining the same concentrations of targets A and B, forward primer A,and reverse primers A and B, and a variable concentration of forwardprimer B for target B, in accordance with aspects of the presentdisclosure.

FIG. 12 is a schematic view of a partition from an exemplary multiplexeddigital assay performed with two detection channels, with the partitioncontaining a copy of a masked target A and a copy of a masking target B,each amplifiable at a different efficiency with a different pair offorward and reverse primers (F_(A) and R_(A) or F_(B) and “R_(B)”) anddetectable via different probes (P_(A) and P_(B)) having differentfluorophore labels, in accordance with aspects of the presentdisclosure.

FIG. 12A is a schematic scatter plot of amplification data that may becollected from the multiplexed digital assay of FIG. 12, with threedistinct clusters of partitions numbered sequentially and identifiedaccording to target content, in accordance with aspects of the presentdisclosure.

FIG. 13 is a schematic view of a partition from another exemplarymultiplexed digital assay performed with two detection channels, withthe partition containing a copy of a masked target (A) and a copy of amasking target (B), each amplifiable with a different pair of forwardand reverse primers (F_(A) and R_(A) or F_(B) and R_(B)) and detectablevia a generic reporter (e.g., an intercalating dye) and a specific probe(P_(B)), in accordance with aspects of the present disclosure.

FIG. 14 is a schematic scatter plot of amplification data that may becollected from the multiplexed digital assay of FIG. 13, with threedistinct clusters of partitions numbered sequentially and identifiedaccording to target content, in accordance with aspects of the presentdisclosure.

FIG. 14A is a schematic view of a partition from an exemplarymultiplexed digital assay performed with a single detection channel,with the partition containing a copy of a masked target (A) and a copyof a masking target (B), each amplifiable with a different pair offorward and reverse primers (F_(A) and R_(A) or F_(B) and R_(B)) anddetectable via a generic reporter (e.g., an intercalating dye) and aspecific probe (P_(B)), in accordance with aspects of the presentdisclosure.

FIG. 14B is a schematic graph of fluorescence intensity data that may becollected from five sets of partitions (lanes 1-5) in a single detectionchannel for the multiplexed digital assay of FIG. 14A, with each setcontaining the same concentrations of targets A and B and of both pairsof primers, and a variable concentration of the probe (P_(B)) for targetB, in accordance with aspects of the present disclosure.

FIG. 15 is a schematic view of a partition from another exemplarymultiplexed digital assay performed with two detection channels, withthe partition containing a copy of a masked target (A), a second target(B), and a masking target (C), each amplifiable with a different pair offorward and reverse primers (F_(A) and R_(A), F_(B) and R_(B), or F_(C)and R_(C)) and detectable via a generic reporter (e.g., an intercalatingdye) and specific probes (P_(B) and P_(C)), in accordance with aspectsof the present disclosure.

FIG. 16 is a schematic scatter plot of amplification data that may becollected from the multiplexed digital assay of FIG. 15 in a firstchannel that detects amplification of targets A and C and in a secondchannel that detects amplification of target B, with five distinctclusters of target-positive partitions numbered sequentially andidentified according to target content, in accordance with aspects ofthe present disclosure.

FIG. 17 is another schematic scatter plot of amplification data that maybe collected from the multiplexed digital assay of FIG. 15 and plottedas in FIG. 16, except that two of the clusters ((A) and (AB)) of FIG. 16are no longer resolved and form a heterogeneous cluster withheterogeneous target content.

FIG. 18 is another schematic scatter plot of amplification data that maybe collected from the multiplexed digital assay of FIG. 15, but with afirst channel that detects amplification of targets A and C and a secondchannel that detection amplification of targets A and B, with fourdistinct clusters of target-positive partitions numbered sequentiallyand identified according to target content, in accordance with aspectsof the present disclosure.

FIG. 19 is a schematic view of a partition from an exemplary multiplexeddigital assay to determine the quality of a library, with the partitioncontaining a copy of an empty library member (an adapter-adapterinverted repeat with no insert) and a copy of a library membercontaining an insert, with each library member being amplifiable withtwo copies of the same primer (F) and detectable with anadapter-specific probe (P), in accordance with aspects of the presentdisclosure.

FIG. 20 is a graph of amplification data collected from the multiplexedassay of FIG. 19 performed in droplets.

FIG. 21 is a schematic view of a partition from an exemplary multiplexeddigital assay to quantify spliced and unspliced species in a sample,with the partition containing a copy of an unspliced species(Exon-Intron-Exon) and a spliced species (Exon-Exon), with each speciesbeing amplifiable with the same pair of forward and reverse primers (Fand R) and detectable via a generic reporter, in accordance with aspectsof the present disclosure.

FIG. 22 is a graph of amplification data collected from the multiplexedassay of FIG. 21 performed in droplets.

FIG. 23 is a schematic view of a partition from another exemplarymultiplexed digital assay to quantify spliced and unspliced species in asample, with the partition containing a copy of an unspliced species(Exon-Intron-X-Intron-Exon) and two distinct spliced species(Exon-X-Exon and Exon-Exon), with each species being amplifiable withthe same pair of forward and reverse primers (F and R) and detectablevia a generic reporter, in accordance with aspects of the presentdisclosure.

FIG. 24 is a graph of amplification data collected from the multiplexedassay of FIG. 23.

FIG. 25 is a schematic view of a partition from an exemplary multiplexeddigital assay to quantify a pair of targets (targets A and B) inpartitions, with the depicted partition containing a copy of each targetand only a single reporter that binds to only one of the targets (targetA), in accordance with aspects of the present disclosure.

FIG. 26 is a schematic plot of amplification data that may be collectedfrom the multiplexed digital assay of FIG. 25 in a pair of opticalchannels, with only one of the channels detecting the single reporter,and with the plot having three distinct clusters of partitions numberedsequentially (0, 1, and 2) and identified according to target content,in accordance with aspects of the present disclosure.

DETAILED DESCRIPTION

The present disclosure provides a multiplexed digital assay system,including methods and apparatus, for calculating a level of one or moretargets. In an exemplary method, data for amplification of a pluralityof targets including a first target may be collected from partitions. Alevel of the first target may be calculated from only a subset of thedata that excludes partitions according to target content.

An exemplary method of performing a multiplexed digital assay isprovided. In the method, a mixture may be provided. The mixture mayinclude a plurality of targets and reagents sufficient for amplificationof each of the targets, with the plurality of targets including a firsttarget. Partitions may be formed, with each partition including aportion of the mixture. The plurality of targets may be amplified in thepartitions. Data for amplification of each of the plurality of targetsmay be collected from a plurality of the partitions. A level of thefirst target may be calculated from only a portion of the data. Theplurality of targets also may include a second target. The portion ofthe data may or may not exclude each partition testing positive for thesecond target.

In some embodiments, partition populations may be identified in thedata, with each partition population having a different target contentof the plurality of targets. The portion of the data may be collectedfrom only a subset of the plurality of the partitions, and the subset ofpartitions may exclude two or more of the partition populations. The twoor more excluded partition populations may include at least a pair ofexcluded partition populations that are not resolved from each other inthe data. The at least a pair of excluded populations may beheterogeneous with respect to the first target.

In some embodiments, the data may be plottable to form three or moreclusters of data points. Each data point may represent a partition. Eachcluster may have a different content of the plurality of targets. Thelevel of the first target may be calculated from only a portion of thedata that excludes at least one of the clusters of data points. The atleast one excluded cluster may include a heterogeneous excluded clustercomposed of overlapping sets of data points representing at least twoexcluded partition populations. At least one of the excluded partitionpopulations may be positive for the first target, and at least one ofthe excluded partition populations may be negative for the first target.

The digital assays of the present disclosure may have numerousadvantages, such as higher levels of multiplexing, multiplexing with ageneric reporter alone or in combination with at least one specificreporter, determination of target levels (e.g., concentrations) withpartition populations that are indistinguishable and/or not wellresolved in the data, determination of target levels using competingassays within partitions, or any combination thereof, among others.

Further aspects of the present disclosure are presented in the followingsections: (I) overview of multiplexed digital assays with a genericreporter, (II) determination of concentration for masked targets, (III)exemplary assay configurations with target masking, (IV) multiplexedassays with a specific reporter and a generic reporter, and (V)examples.

I. OVERVIEW OF MULTIPLEXED DIGITAL ASSAYS WITH A GENERIC REPORTER

This section provides an overview of multiplexed digital assaysperformed with a generic reporter and, optionally, in which the presenceof a target in a partition is at least partially masked by the presenceof a different target in the same partition; see FIGS. 1 and 2. Thedigital assays may utilize data exclusion as described Section II andelsewhere herein to calculate a level of one or more of the targets.

FIG. 1 shows a flowchart of an exemplary method 50 of performing amultiplexed digital assay with a generic reporter, optionally withtarget masking. The steps presented for method 50 may be performed inany suitable order and in any suitable combination. Furthermore, thesteps may be combined with and/or modified by any other suitable steps,aspects, and/features of the present disclosure, including thosedescribed in the patent documents listed above under Cross-References,which are incorporated herein by reference.

Sample Preparation.

A sample may be prepared for the assay, indicated at 52. Preparation ofthe sample may include any suitable manipulation of the sample, such ascollection, dilution, concentration, purification, lyophilization,freezing, extraction, combination with one or more assay reagents toform a mixture (also termed a sample-containing mixture, a bulk phase,or a reaction mixture), performance of at least one preliminary reactionto prepare the sample for one or more reactions in the assay, or anycombination thereof, among others. The preparation may isolate ananalyte, such as nucleic acid that includes copies of one or morenucleic acid targets, and/or may modify and/or fragment the analyte.Preparation of the sample may include rendering the sample competent forsubsequent performance of one or more reactions, such as one or moreenzyme catalyzed reactions and/or binding reactions.

In some embodiments, preparation of the sample may include combining thesample with reagents to produce a sample-containing mixture forperforming a reaction (such as an amplification reaction) for eachtarget and for reporting an extent of each reaction (e.g., whether ornot the reaction occurred above a threshold level or within a range).Reagents for amplification may include any combination of primers fortargets, dNTPs and/or NTPs, at least one enzyme (e.g., a polymerase, aligase, a reverse transcriptase, a restriction enzyme, or a combinationthereof, among others, each of which may or may not be heat-stable),and/or the like. Accordingly, the mixture may have a complete set ofreagents for (i.e., may be competent for) amplification of each targetunder suitable environmental conditions (e.g., incubation at an elevatedtemperature or modulation of temperature (such as by thermocycling)).The mixture may be capable of amplification of each of one or moretargets, if present, in the sample (or a partition thereof). Reagentsfor reporting may include at least one generic reporter and/or at leastone specific reporter. The generic reporter may be sensitive toamplification of each target and the specific reporter may bespecifically sensitive to amplification of only a subset of the targets,such as only one of the targets. The mixture may or may not include adifferent reporter for each target to be assayed. Preparation of themixture may render the sample capable of reporting, or being analyzedfor, whether or not a reaction, such as amplification, has occurred, ona target-by-target basis, and optionally the extent of any suchreaction.

Providing Partitions.

Partitions for the assay may be provided, indicated at 54. Eachpartition may include a portion of a same mixture. In some cases, theportion may constitute the entire partition. The mixture may containeach target (e.g., provided by a same sample), each reporter, and/or oneor more amplification reagents (e.g., a complete set of reagents foramplification of each target). Accordingly, the partitions,collectively, may contain a plurality of targets and each partition maycontain the same generic reporter and optionally the same specificreporter. The targets may include at least one maskable target(interchangeably termed a masked target), and at least one obscuringtarget (interchangeably termed a masking target or a dominant target).The obscuring target is capable of at least partially masking thepresence of one or maskable targets when a copy of each target ispresent in the same partition.

The partitions when provided (e.g., when formed) may contain each targetat “partial occupancy,” which means that each partition of only a subsetof the partitions contains at least one copy of each target (and/ortemplate) to be assayed. For example, with a multiplexed assay performedon a first target and a second target, only a first subset of thepartitions contains the first target, and only a second subset of thepartitions contains the second target. The first subset and the secondsubset of the partitions may be the same subset, if the first target andthe second target are fully associated with and/or linked to each otherwhen the partitions are formed. Alternatively, the first subset and thesecond subset of the partitions may be distinct if the first target andthe second target are not fully associated with and/or linked to eachother when the partitions are formed. In some cases, if the targets arenot fully associated and/or linked, each partition of a distinct thirdsubset of the partitions may contain at least one copy of each target.Accordingly, with partial occupancy, one or more (e.g., a plurality) ofthe partitions contain no copies of the first target, one or more (e.g.,a plurality) of the partitions may contain a single copy (only one copy)of the first target, and, optionally, yet one or more of the partitions(e.g., the rest of the partitions) may contain two or more copies of thefirst target. Similarly, with partial occupancy, one or more (e.g., aplurality) of the partitions contain no copies of the second target, oneor more (e.g., a plurality) of the partitions may contain a single copy(only one copy) of the second target, and, optionally, yet one or moreof the partitions (e.g., the rest of the partitions) may contain two ormore copies of the second target.

The term “partial occupancy” is not restricted to the case where thereis no more than one copy of a particular template/target of interest inany partition. Partitions containing a template and/or a target atpartial occupancy may, for example, contain an average of more than, orless than, about one copy, two copies, or three copies, among others, ofthe template/target per partition when the partitions are provided orformed. Copies of a template (and/or target) may have a randomdistribution among the partitions, which may be described as a Poissondistribution. In some cases, a significant number of the partitions(e.g., at least about 1%, 2%, 5%, 10%, or 20%, among others, of thepartitions) may contain a copy of each of at least two targets, and/or aplurality of the partitions each may contain at least one copy of alltargets.

Targets may be unlinked, partially linked, or fully linked whenpartitions are formed. Linked targets may be attached to each othercovalently and/or by base pairing.

Each target may be an analyte on interest (e.g., a nucleotide sequenceof interest) or a surrogate therefor (e.g., a nucleic acid bound toand/or corresponding to a nucleotide sequence of interest). The targetmay be nucleic acid that includes a sequence of nucleotides. The target,if nucleic acid, may be single-stranded or double-stranded, amongothers. A nucleic acid target may be provided by a template, with thetarget forming at least a portion or all of the template. The target maycorrespond to an amplicon produced by amplification. The amplicon may besingle-stranded or double-stranded, among others. In some cases, thetarget may be or correspond to an analyte that is not nucleic acid, suchas a small molecule, a polypeptide, a lipid, an amino acid, an ion, orthe like.

The partitions may be provided by distributing or separating portions ofa sample-containing bulk phase into partitions. Any suitable fractionincluding up to all of the bulk phase may be distributed to thepartitions. Each partition may be and/or include a fluid volume that isisolated from the fluid volumes of other partitions. The partitions maybe isolated from one another by a fluid phase, such as a continuousphase of an emulsion, by a solid phase, such as at least one wall of acontainer, or a combination thereof, among others. In some embodiments,the partitions may be droplets disposed in a continuous phase, such thatthe droplets and the continuous phase collectively form an emulsion.

The partitions may be formed by any suitable procedure, in any suitablemanner, and with any suitable properties. For example, the partitionsmay be formed with a fluid dispenser, such as a pipette, with at leastone droplet generator, by agitation of the sample (e.g., shaking,stirring, sonication, etc.), and/or the like. Accordingly, thepartitions may be formed serially, in parallel, or in batch. Thepartitions may have any suitable volume or volumes. The partitions maybe of substantially uniform volume or may have different volumes.Exemplary partitions having substantially the same volume aremonodisperse droplets. Exemplary volumes for the partitions include anaverage volume of less than about 100, 10 or 1 μL, less than about 100,10, or 1 mL, or less than about 100, 10, or 1 μL, among others.

The partitions, when initially formed, may be competent for performanceof one or more reactions in the partitions. Alternatively, one or morereagents may be added to the partitions after they are formed, toprovide partitions competent for reaction. The reagents may be added byany suitable mechanism, such as a fluid dispenser, fusion ofpartitions/droplets, or the like.

Performance of Reactions.

A reaction corresponding to each target may be performed in thepartitions, indicated at 56. The reaction may be an enzyme-catalyzedreaction. The reaction may produce a product, which may increaselinearly or exponentially, among others, during the reaction.

The reaction for each target may be an amplification reaction thatamplifies copies of the target present in individual partitions.Amplification of each target may occur selectively in partitionscontaining at least one copy of the target (e.g., containing at leastone copy of a template that includes the target). Amplification may belinear or exponential, among others.

Amplification may or may not be performed isothermally. In some cases,amplification in the partitions may be encouraged by heating thepartitions and/or incubating the partitions at a temperature above roomtemperature, such as at a denaturation temperature, an annealingtemperature, and/or an extension temperature, for one or a plurality ofcycles. In some examples, the partitions may be thermally cycled topromote a polymerase chain reaction and/or ligase chain reaction.Exemplary isothermal amplification approaches that may be suitableinclude nucleic acid sequence-based amplification,transcription-mediated amplification, multiple-displacementamplification, strand-displacement amplification, rolling-circleamplification, loop-mediated amplification of DNA, helicase-dependentamplification, and single-primer amplification, among others.

Data Collection.

Data may be collected from the reporters in the partitions, indicated at58. The partitions may contain only one reporter or a plurality ofdifferent reporters having any suitable structure and characteristics.Each reporter may be a specific reporter (e.g., a probe) or a genericreporter (e.g., an intercalating dye). In some cases, at least onespecific reporter may report on the occurrence of only one particularreaction and thus only one target in a multiplexed assay. In othercases, a specific reporter may report on the occurrence of two or morereactions and thus two or more targets. Each reporter interchangeablymay be termed a reaction reporter and/or an amplification reporter.

A specific reporter binds with substantial specificity to a product(e.g., an amplicon) of a reaction, based on an identity of the product,to the substantial exclusion of other structurally different substancesof the same class as the product. For example, each specific reportermay bind to a particular sequence or site to the substantial exclusionof other sequences or sites. The specific reporter may include or be aprobe including at least one nucleic acid (e.g., at least oneoligonucleotide) that binds specifically to a complementary nucleotidesequence present in a target and/or an amplicon corresponding to thetarget. The probe may include a label associated with an oligonucleotide(e.g., covalently attached or bound noncovalently, among others). Thelabel may be optically detectable directly or indirectly. Accordingly,the label may be a luminophore (such as a photoluminescent moiety (e.g.,a fluorophore or phosphor)), an enzyme (e.g., a peroxidase,beta-galactosidase, alkaline phosphatase, phosphodiesterase, or thelike), a member of a specific binding pair (e.g., biotin oravidin/streptavidin), or an epitope tag, among others. A probe includinga luminophore may or may not also include an energy transfer partner forthe luminophore, such as a quencher or another luminophore (e.g., toproduce luminescence resonance energy transfer (e.g., FRET)). The probemay or may not also function as a primer that is extended in the assay.Exemplary labeled probes include TaqMan® probes, Scorpion®probes/primers, Eclipse® probes, Amplifluor® probes, molecular beaconprobes, Lux® primers, proximity-dependent pairs of hybridization probesthat exhibit FRET when bound adjacent one another on an amplicon, QZyme®primers, or the like.

The specific reporter may have distinct forms or states. The specificreporter may have an initial/intact form or state and one or moredegraded/modified forms or states. The one or more degraded/modifiedforms or states may be produced from the initial/intact form duringamplification of the second target. The forms or states may bedistinguishable optically. For example, the degraded/modified form maybe more or less photoluminescent than the initial/intact form.

A generic reporter (interchangeably termed a nonspecific reporter) bindswithout substantial specificity to a product of a reaction (e.g., anamplicon), such that other structurally different substances of the sameclass as the product (e.g., other amplicons of unrelated sequence) canalso be bound by the reporter. The nonspecific binding may not depend ona unique feature of the arrangement of atoms of one or both of thereporter and the product (e.g., the target and/or amplicon). Multiplecopies of the generic reporter may be capable of binding to a singlecopy of a reaction product, for example, with the number of copies boundbeing related directly, such as proportional, to the amount or length ofthe reaction product. For example, the generic reporter may be aphotoluminescent dye that binds to nucleic acid relativelynonspecifically. The dye may not be attached to an oligonucleotide thatconfers substantial sequence binding specificity. The dye may be a majorgroove binder, a minor groove binder, an intercalator, or an externalbinder, among others. The dye may bind preferentially to double-strandedrelative to single-stranded nucleic acid and/or may exhibit a greaterchange in a photoluminescence characteristic (e.g., emission intensity)when bound to double-stranded relative to single-stranded nucleic acid.Exemplary dyes that may be suitable include luminescent cyanines,phenanthridines, acridines, indoles, imidazoles, and the like, such asDAPI, Hoechst® 33258 dye, acridine orange, etc. Exemplary intercalatingdyes that may be suitable include ethidium bromide, propidium iodide,EvaGreen® dye, SYBR® Green dye, SYBR® Gold dye, and 7-aminoactinomycin D(7-AAD), among others.

R targets may be assayed in a multiplex assay, and the data may becollected in less than R optical channels (e.g., in differentwavebands). In other words, the number (R) of targets assayed may begreater than the number of optical channels used for detecting thetarget-specific reactions. In some cases, the data may be collected inonly one or two optical channels, or in at least two, three, or moreoptical channels, among others. In some cases, data may be collectedfrom the same number of optical channels as targets in the assay. Anoptical channel interchangeably may be termed a detection channel.

An optical channel may represent a particular detection regime withwhich emitted light is generated and detected. The detection regime maybe characterized by a spectral content (i.e., a wavelength regime) fordetection of emitted light. If pulsed excitation light is used in thedetection regime to induce light emission, the detection regime may becharacterized by a spectral content (a wavelength(s) or waveband(s)) forillumination with excitation light and/or a time interval during whichlight emission is detected with respect to each light pulse.Accordingly, optical channels that are different from each other maydiffer with respect to the spectral content (wavelength(s)/waveband(s))of excitation light, with respect to the spectral content(wavelength(s)/waveband(s)) of emitted light that is detected, and/orwith respect to the time interval during which emitted light is detectedrelative to each pulse of excitation light, among others.

Data collection may include generating one or more signalsrepresentative of detected light. The signal may represent an aspect oflight, such as the intensity of the light, detected in the same opticalchannel from one or more reporters for two or more distinct targets. Thesignals optionally may include data collected in two or more differentoptical channels (e.g., at different wavelengths and/or differentwavelength ranges (wavebands) and/or color regimes) from reporters forthe same and/or different targets). The light detected from eachreporter may be light emitted from a luminophore (e.g., a fluorophore).The light detected in a given channel may be detected such that lightfrom different reporters is summed or accumulated without attribution toa particular reporter. Thus, the signal for a given channel may be acomposite signal that represents two, three, four, or more reportersand/or assays and two, three, four, or more targets.

The signal(s) may be created based on detected light emitted by thegeneric reporter, and optionally from one or more other reporters in thepartitions. The one or more other reporters may report whether at leastone of two or more particular amplification reactions represented by thesignal has occurred in a partition and thus whether at least one copy ofat least one of two or more particular targets corresponding to the twoor more particular amplification reactions is present in the partition.The level or amplitude of the signal corresponding to the reporters maybe analyzed to determine whether or not at least one of the particularreactions has occurred and at least one copy of one of the particulartargets is present. The level or amplitude of the signal may vary amongthe partitions according to whether at least one of the particularreactions occurred or did not occur and at least one of the particulartargets is present or absent in each partition. For example, a partitionpositive for a particular target may produce a signal level or amplitudethat is above a given threshold and/or within a given range. Partitionsmay be analyzed and signals created at any suitable time(s). Exemplarytimes include at the end of an assay (endpoint assay), when reactionshave run to completion and the data no longer are changing, or at someearlier time, as long as the data are sufficiently and reliablyseparated.

Data may be collected from a plurality of the partitions (i.e., only asubset or all of the partitions) under any suitable conditions. All ofthe data may be collected at about the same temperature from theplurality of partitions, at a temperature that is below a meltingtemperature of each amplicon, and/or below about 50 degrees Celsius,among others. The amplification data may be collected after an endpointof amplification has been reached for each target.

Partition clusters each positive for at least one different combinationof zero, one, or more of the targets may be identified. Generally, ifmasking occurs, only a subset of all possible combinations may beidentified and some of the clusters may contain two or more populationshaving distinct target compositions. Identification may be performed bya data processor using an algorithm (e.g., an algorithm that identifiespatterns (e.g., partition clusters) in the data), by a user, or acombination thereof. In some cases, a data processor may produce andoutput (e.g., display) a plot of the collected data (e.g., a 2-D scatterplot or histogram, or, with three or more optical channels fordetection, two or more 2-D scatter plots or histograms with differentpairs of axes). The user then may define the boundary of each clusterbased on the plot(s), e.g., through a graphical user interface to definepopulation boundaries, and/or by inputting values (e.g., representingintensity ranges) to define a boundary for each cluster. Each clusterboundary may be defined by one or more ranges of values, a geometricalshape that surrounds the cluster (e.g., a polygon, ellipse, etc.), orthe like. Further aspects of cluster identification are presented belowin Sections II-V.

The reaction components and/or conditions of any of the multiplexedassays disclosed herein may be adjusted to improve the resolution ofdifferent partition populations in the data. By changing theconcentration of a particular assay within a multiplexed assay, thereaction efficiency for a particular target can be affected, which mayresult in a difference in signal level that allows populations detectedwith the same reporter and/or different reporters to be distinguishedfrom one another. By changing reaction components/conditions, additionaltargets may be detected in the same multiplexed reaction. In some cases,the signal amplitude for a target may be adjusted by varying theconcentration of one or both primers for the target. Varying primerconcentration without changing the reporter concentration may be usefulin assays where the same reporter (e.g., a probe or a generic reporter)is used to detect two or more targets, but each of the two targets isamplified with at least one different primer. In some cases, the signalamplitude for one or more targets may be adjusted by changing theannealing temperature used for thermocycling, the total concentration ofdNTPs, the amounts of individual dNTPs relative to each other (e.g., ifthe two targets have substantially different base compositions), or thelike.

A partition count for each partition cluster may be obtained. Thepartition count may be a value representing the number of partitionsconstituting a particular partition cluster.

A number of partitions that are positive (or negative) for each targetmay be determined from the collected data. The signal detected from eachpartition, and the partition itself, may be classified as being positiveor negative for each of the reactions/targets contributing to thesignal. Classification may be based on the strength (and/or othersuitable aspect) of the signal. For example, classification may be basedon the intensity of light detected from the partition in a singleoptical channel or in two or more optical channels. If thesignal/partition is classified as (testing) positive (+), for a giventarget, the reaction corresponding to that target is deemed to haveoccurred and at least one copy of the target is deemed to be present inthe partition. In contrast, if the signal/partition is classified asnegative (−), for a given target, the reaction corresponding to thattarget is deemed not to have occurred and no copy of the target isdeemed to be present in the partition (i.e., the target is deemed to beabsent from the partition).

The data including all combinations of targets (i.e., all combinationsof the presence or absence of each of the targets) will generallyrepresent 2^(R) populations each having a different target content,where R is the number of different targets. For example, with threetargets (A-C, R=3), there are eight populations: ( ), (A), (B), (C),(AB), (AC), (BC), and (ABC), each having a different target content. Iftarget masking occurs, some of the populations will overlap and/or willnot be distinguishable from one another as well-resolved clusters.Instead, some of the populations may overlap one another to form atleast one heterogeneous cluster having a heterogeneous target content ofone or more of the targets. Accordingly, less than 2^(R) well-resolvedclusters may be produced with target masking and/or with other assayconfigurations as described elsewhere in the present disclosure.

Determination of Target Concentrations.

A concentration of each target may be determined, indicated at 60. Theconcentration of at least one of the targets may be determined from onlya subset of the collected data. For example, the concentration of atleast one masked target may be determined from only a subset of thedata. The subset of data may selectively exclude partitions that arepositive for an obscuring target that at least partially obscures thepresence of the masked target in partitions containing a copy of bothtargets.

Any suitable subset of the partitions may be excluded, for calculatingthe concentration of a given target, without skewing the concentrationdetermined, if the basis for exclusion is independent of thepresence/absence of the target. For example, in an assay of targets Aand B, all (or any suitable subset of) B-positive partitions may beexcluded from calculation of the concentration for target A, if beingpositive for B is independent of the presence or absence or A. (Forexample, there is no substantial association or linkage between targetsA and B when partitions are formed.) Accordingly, if all B-positivepartitions are present in the same cluster (or two or more clusters),the entire cluster (or the two or more clusters) can be excluded fromthe calculation of the concentration of target A, without skewing theresult. As another example, in an assay for unlinked targets A, B, andC, the concentration of target A can be determined with all B-positivesand all C-positives excluded, with only B-positives or only C-positivesexcluded, or with any combination of B-positives and C-positivesexcluded that is independent of the presence/absence of A (e.g., allBC-positives, whether A-negative or A-positive). In the same assay, theconcentration of target B (or target C) may be determined from all ofthe data, if all B-positives (or all C-positives) are resolved from allother populations lacking a copy of target B (or target C).Alternatively, the concentration of target B (or target C) may bedetermined from only a subset of the data that excludes partitionsindependent of the presence/absence of B (or C).

Determination of target concentrations may (or may not) be based on eachtarget having a Poisson distribution among the partitions. Eachconcentration may, for example, be a value representing the averagenumber copies of the target per partition. The partition data furthermay be used (e.g., directly and/or as concentration data) to estimatecopy number (CN) and copy number variation (CNV), using any suitablealgorithms.

A level, such as a concentration, of each target may be determined withPoisson statistics. The concentration of the target in the partitionsmay be calculated from the fraction of partitions that are positive forthe target (or, equivalently, the fraction of partitions that arenegative for the target) by assuming that copies of the target (beforeamplification) have a Poisson distribution among the partitions. Withthis assumption, the fraction f(k) of partitions having k copies of thetarget is given by the following equation:

$\begin{matrix}{{f(k)} = {\frac{\lambda^{k}}{k!}e^{- \lambda}}} & (1)\end{matrix}$

Here, A is the concentration of the target in the partitions, expressedas the average number of target copies per partition (beforeamplification). Simplified Poisson equations may be derived from themore general equation above and may be used to determine targetconcentration from the fraction of positive partitions. An exemplaryPoisson equation that may be used is as follows:

$\begin{matrix}{\lambda = {- {\ln \left( {1 - \frac{N_{+}}{N_{tot}}} \right)}}} & (2)\end{matrix}$

where N₊ is the number of partitions (i.e., the partition count)positive for a given target, and where N_(tot) is the total number ofpartitions that are positive or negative for the target. N_(tot) isequal to a sum of (a) N₊ for the target and (b) the number of partitionsnegative for the target, or N⁻. N₃₀/N_(tot) (or N₊/(N₊+N⁻) is equal tof₊, which is the fraction of partitions positive for the target (i.e.,f₊=f(1)+f(2)+f(3)+ . . . ) (see Equation 1), and which is a measuredestimate of the probability of a partition having at least one copy ofthe target. Another exemplary Poisson equation that may be used is asfollows:

$\begin{matrix}{\lambda = {- {\ln \left( \frac{N_{-}}{N_{tot}} \right)}}} & (3)\end{matrix}$

where N⁻ and N_(tot) are as defined above. N⁻/N_(tot) is equal to f₀,which is the fraction of negative partitions (or 1−f₊), is a measuredestimate of the probability of a partition having no copies of thetarget, and A is the target concentration as described above.

Equations 2 and 3 above can be rearranged to produce the following:

λ=ln(N _(tot))−ln(N _(tot) −N ₊)  (4)

Equations 2 and 3 above can be rearranged to produce the following:

λ=ln(N _(tot))−ln(N ⁻)  (5)

The concentration of each target in a multiplexed assay can, forexample, be determined with any of Equations 2-5, using values (i.e.,partition counts) obtained for N_(tot) and N⁻ or N₊, for each target. Insome cases, the value used for N_(tot) (the total partition count) maybe the same for each target. In other cases, the value used for N_(tot)may vary, such as if some of the populations are excluded from the totalcount due to population overlap (target masking). In some embodiments,N_(tot) may be equivalent to a combination of all populations, namely, asum of the partition counts for all populations identified.

The value used for N⁻ or N₊ is generally different for each target, andmay result from summing the counts from a plurality of partitionpopulations each containing a different combination of the targets beingtested in the multiplexed assay. For example, with three targets (A, B,and C) in a multiplexed assay, the number of partitions positive fortarget A, N_(+A), may be calculated as the sum of counts from the single(A only), double (AB and AC), and triple (ABC) positive populations, foruse in Equation 2 or 4. Equivalently, the number of partitions negativefor target A, N_(−A), may be calculated, for use in Equation 3 or 5, asthe difference between N_(tot) and N_(+A). Alternatively, the number ofpartitions negative for A may be calculated as the sum of counts fromeach population that is negative for target A, namely, in this example,a triple negative (“empty”) population, two single positive populations(B and C), and one double positive population (BC). The same process maybe repeated for each of the other targets using partition counts fromthe appropriate subset of populations.

In some embodiments, an estimate of the level of the target may beobtained directly from the positive fraction, without use of Poissonstatistics. In particular, the positive fraction and the concentration(copies per partition) converge as the concentration decreases. Forexample, with a positive fraction of 0.1, the concentration isdetermined with Equation 2 to be about 0.105, a difference of only 5%;with a positive fraction of 0.01, the concentration is determined to beabout 0.01005, a ten-fold smaller difference of only 0.5%. However, theuse of Poisson statistics can provide a more accurate estimate ofconcentration, particularly with a relatively higher positive fraction,because the equation accounts for the occurrence of multiple targetcopies per partition.

FIG. 2 shows an exemplary system 70 for performing the digital assay ofFIG. 1. System 70 may include a partitioning assembly, such as a dropletgenerator 72 (“DG”), a thermal incubation assembly, such as athermocycler 74 (“TC”), a detection assembly (a detector) 76 (“DET”),and a data processing assembly (a data processor) 78 (“PROC”), or anycombination thereof, among others. The data processing assembly may be,or may be included in, a controller that communicates with and controlsoperation of any suitable combination of the assemblies. The arrowsbetween the assemblies indicate movement or transfer of material, suchas fluid (e.g., a continuous phase of an emulsion) and/or partitions(e.g., droplets) or signals/data, between the assemblies. Any suitablecombination of the assemblies may be operatively connected to oneanother, and/or one or more of the assemblies may be unconnected to theother assemblies, such that, for example, material/data are transferredmanually.

Detector 76 may provide a plurality of optical channels in which datacan be collected. The detector may have a distinct sensor or detectionunit for each optical channel.

System 70 may operate as follows. Droplet generator 72 may form dropletsdisposed in a continuous phase. The droplets may be cycled thermallywith thermocycler 74 to promote amplification of targets in thedroplets. Signals may be detected from the droplets with detector 76.The signals may be processed by processor 78 to determine numbers ofdroplets and/or target concentrations, among others. The system mayinclude a program, optionally residing on a computer-readable storagemedium, and comprising instructions for causing the data processorand/or controller to perform and/or control any suitable combination ofthe steps disclosed herein, such as in FIG. 1.

Further aspects of sample preparation, partition formation (such asdroplet generation), data collection, population identification and/orcluster assignment, obtaining partition counts, and target leveldetermination, among others, that may be suitable for the system of thepresent disclosure are described in Sections II-V and in the referencesidentified above under Cross-References, which are incorporated hereinby reference.

II. DETERMINATION OF CONCENTRATION FOR MASKED AND MASKING TARGETS

This section describes exemplary approaches to determining targetconcentrations when target masking occurs and/or when partitionpopulations overlap with one another in the collected data; see FIGS.3-5.

In a typical duplex (two-target) assay performed with two or moreoptical channels for detection, the presence of different targets in thesame partitions leads to signals in different channels. However, in somecases, only one detection channel may be utilized a duplex assay, whichmay cause one target to mask the detectability of another. For example,masking can occur when two amplicons of significantly differentamplification efficiencies are produced in the same droplet and aredetectable via the same generic reporter, such as a SYBR®-derived dye.The presence of the amplicon that accumulates more slowly can be maskedby the amplicon that accumulates more quickly (e.g., see Example 2).Another example can be found in next-generation sequence (NGS) libraryformation. Products that only contain a pair of adapters without aninsert may outcompete fully formed species and thus mask them in asingle-channel assay (e.g., see Example 1). The present disclosureenables robust calculation of the concentration of a pair of competingtargets that produce target masking.

In a single-channel assay, it may be possible to accurately measure theconcentrations of multiple competing species provided: (1) the speciesproduce different signal levels, (2) rank order masking or suppressionof species exists, and (3) the signal is separated into a bandingpattern allowing counting of partitions in bands. Examples of suchsituations in practice include variable size amplicons, such as mRNAversus precursor RNA, loci with small insertions or deletions, NGSlibrary quantitation, or the like.

FIG. 3 shows a schematic graph of amplification data that may becollected in a single detection channel. The graph plots signalamplitude (e.g., photoluminescence intensity, such as fluorescenceintensity) as a function of detected event number. Each event representsa spike or wave in signal strength caused by a detected partition 90(such as a droplet 92 flowing through a detection region). The amplitudeplotted for each partition/event may be a maximum signal value, anaverage signal value, or an integrated signal value, among others, forthe corresponding spike or wave. The x-axis of FIG. 3 alternatively maybe labeled as time or partition number, among others.

The partitions represented by the data of FIG. 3 collectively contain amasking target 94 (“H”) and a masked target 96 (“L”) that producerespective, characteristic higher and lower amplitude signals. When atleast one copy of target H is present in a partition, the higheramplitude signal is dominant, as exemplified at 98 by a group ofH-positive partitions. The higher amplitude signal is detected from eachH-positive partition, whether or not a copy of target L also is presentin the partition. In other words, the signal produced by amplificationof target H masks the signal from amplification of target L, becausetarget L does not sufficiently increase (or decrease) the signalamplitude over that generated by amplification of target H alone. Incontrast, when one or more copies of target L alone are present in apartition (i.e., in the absence of any copies of target H), the loweramplitude signal is detected from the partition, as exemplified at 100by a group of L-positive, H-negative partitions. Individual partitions,such as partition 102 that do not contain a copy of either target(“empties” or “negatives”; abbreviated “E”) produce an even lower,baseline signal. Accordingly, three bands 104, 106, 108 of partitionsare visible in the graph. (Bands interchangeably may be termedclusters.) Band 104 represents H-positive partitions (±L), band 106represents L-positive partitions that are H-negative, and band 108represents H- and L-negative partitions.

The number of partitions in each band can be determined (interchangeablytermed counted), and the concentration of each target calculated using,for example, Equation 5 of Section I. However, the concentration ofmasked target L will be underestimated with this approach, if all of thepartitions are included in the calculation. Accordingly, only a subsetof the data can be utilized to determine the concentration of target L.

The general formula for calculating a concentration value (copies perpartition (cpp)) of a particular factor (e.g., a single species ortarget, a combination of species or targets, etc.) is as follows:

λ_(factor)=ln(N _(tot))−ln(N _(not factor))  (6)

where λ_(factor) is the concentration in cpp of the factor, N_(tot) isthe total number of partitions observed, and N_(not factor) is thenumber of partitions that do not exhibit the factor (which is written asN⁻ in Equations 3 and 5).

Since all H-positive partitions are detectable, the following equationholds:

λ_(H)=ln(N _(E) +N _(L) +N _(H))−ln(N _(E) +N _(L))  (7)

where λ_(H) is the concentration of the masking target H, N_(E) is thenumber of empty partitions, N_(H) is the number of higher amplitudepartitions, and N_(L) is the number of lower amplitude partitions.

The presence of masked target L is independent of the presence ofmasking target H in a partition. Accordingly, the H-positive partitionsare unbiased with respect to target L. In other words, the H-positivepartitions can be excluded from the concentration calculation for maskedtarget L. Stated differently, the concentration for target L computedfrom lower bands 106, 108 should be the same, absent statisticalvariation, as the concentration computed from all three bands (if the Ltarget were not masked). Another way to explain this approach is that ifthe higher amplitude partitions for some reason became invisible, andonly the lower amplitude and empty (negative) bands were visible, theconcentration of target L still could be computed accurately as follows:

λ_(L)=ln(N _(E) +N _(L))−ln(N _(E))  (8)

The total target concentration (i.e., the concentration of L or H) canbe computed in two different ways. A first approach is as follows:

λ_(L or H)=ln(N _(E) +N _(L) +N _(H))−ln(N _(E))  (9)

Alternatively, the two concentrations of targets L and H, calculatedseparately, can be summed:

λ_(L or H)=λ_(H)+λ_(L)  (10)

Substitution with terms from Equations 7 and 8 produces the following:

λ_(L or H)=[ln(N _(E) +N _(L) +N _(H))−ln(N _(E) +N _(L))]+[ln(N _(E) +N_(L))−ln(N _(E))]  (11)

Equation 11 simplifies to the following:

λ_(L or H)=ln(N _(E) +N _(L) +N _(H))−ln(N _(E))  (12)

Therefore, both approaches give the same answer. In any event, moregenerally, the concentration of each masked target can be calculatedfrom only a subset of the partitions that selectively excludes eachpopulation of partitions positive for an obscuring target that at leastpartially obscures the presence of the masked target.

FIG. 3A shows another schematic graph of data that may be collected asin FIG. 3. Here, a population 112 of double-positive partitions (H+L)exhibits a small, detectable increase in signal amplitude over apopulation 114 of single-positive partitions containing only target H.However, the increase in signal amplitude produced by a copy of the Ltarget in population 112 is insufficient to fully resolve the two higheramplitude populations from each other. In other words, the partitions ofeach type, H+L and H-only, cannot be reliably distinguished from eachother; the presence of target H in a partition at least partially masksthe presence of target L. Accordingly, the concentration of target L canbe determined more accurately by excluding both populations 112 and 114from the set of partitions used for the calculation, as in Equation 8above.

FIG. 4 shows another schematic graph of data that may be collectedgenerally as in FIG. 3. Here, however, the masking target is target L,which has a lower amplitude signal, and the masked target is target H,which has a higher amplitude signal. As a result, detectably H-positivepartitions 120 do not contain a copy of target L, while L-positivepartitions 122 can contain no copies, or one or more copies of target H,in the same band or cluster of partitions.

The concentration of masked target H in FIG. 4 can be calculated in amanner analogous to Equation 8 above by excluding the L-positivepartitions from the partition counts used in the calculation:

λ_(H)=ln(N _(E) +N _(H))−ln(N _(E))  (13)

The concentration of target L in FIG. 4 can be calculated in a manneranalogous to Equation 7 above, since target L is not masked and allL-positive partitions are detectable:

λ_(L)=ln(N _(E) +N _(L) +N _(H))−ln(N _(E) +N _(H))  (14)

FIG. 5 shows a schematic graph of data that may be collected in a singledetection channel from partitions, such as droplets, collectivelycontaining k targets that mask one another serially according to amasking hierarchy. More particularly, the k targets can be rank-orderedfrom 1 to k such that the presence of a target of rank i in a partitionmasks the presence of each target with a lower rank (i.e., i-1, i-2, . .. , 1). Here, the signal from a target of rank i is stronger than fromeach target of lower rank, but any assignable order of signal amplitudesmay be suitable. For example, in some cases, the signal amplitude ofeach successive target of higher rank may decrease relative to eachtarget of lower rank (such that target k has a signal amplitude closestto the empty partitions), or the signal amplitudes may not monotonicallyincrease or decrease according to the rank order.

The data for the partitions is distributed among k target-positive bands(B₁, B₂, . . . , B_(k)) plus a band of empty/negative partitions oflowest amplitude (B₀). For each target-positive band of rank i, onlypartitions containing targets with rank i or lower, namely, i-1, i-2, .. . are present in the band (otherwise such a partition would move to aband of higher rank).

The correct concentration of all targets can be calculated using aninductive method. The concentration of each target can be “peeled off”with the following expressions, wherein N₁ is the number of partitionsin band B₁:

λ₁=ln(N ₀ +N ₁)−ln(N ₀)  (15)

λ₂=ln(N ₀ +N ₁ +N ₂)−ln(N ₀ +N ₁)  (16)

And for the i^(th) case:

λ_(i)=ln(N ₀ +N ₁ + . . . +N _(i))−ln(N ₀ +N ₁ + . . . +N _(i-1))  (17)

This approach can be verified as in Equations 9 to 12 above:

λ₁+λ₂+ . . . +λ_(k)=ln(N ₀ +N ₁ + . . . N _(k))−ln(N ₀)  (18)

λ_(tot)=ln(N _(tot))−ln(N ₀)  (19)

The single channel examples presented here may be extended to datacollected in two or more detection channels (e.g., by detection of twoor more distinct wavebands of emitted light in the channels; seeSections III and IV). In some cases, if the data are collected fromchannel-specific reporters, the data may be “collapsed” into a singlechannel for two or more targets that are multiplexed in that channel, byignoring distinguishable populations produced by amplitude differencesin the other detection channel(s).

III. EXEMPLARY ASSAY CONFIGURATIONS WITH TARGET MASKING

This section describes exemplary multiplexed assay configurations inwhich target masking can occur and/or partition populations of differenttarget content in the data; see FIGS. 6-17. Data from any of themultiplexed assays can be processed as described in Section II, usingdata exclusion for calculating at least one level of a target.

Target masking may occur by any suitable mechanism. In some cases,masking may result from different efficiencies of amplification for two(or more) distinct targets in the same partition, with the moreefficiently-amplified target (the dominant target) masking the lessefficiently-amplified target (the maskable target). Competition betweenassays for dominant and maskable targets in the same partition mayattenuate or eliminate the signal change normally detected for assay ofthe maskable target in a partition, in the absence of the dominanttarget. In other words, the competition may prevent the maskable targetassay from reaching its usual endpoint. The dominant target, byoutcompeting the maskable target, may deplete or monopolize a limitingreagent required for the assay of both targets. The limiting reagent maybe a primer used for amplification of both targets, a reporter thatbinds to amplicons for both targets, one or more dNTPs (and/or NTPs), anamplification enzyme that loses substantial activity during the assay,or the like. Alternatively, or in addition, the dominant target, byoutcompeting the maskable target, may produce an inhibitory product,such as nucleic acid (e.g., an amplicon) at a sufficient concentrationto inhibit amplification, to prevent the maskable target assay fromreaching its usual endpoint. In some cases, masking may occur withoutsubstantial competition between targets.

FIG. 6 shows a partition 140 from an exemplary multiplexed digital assaycapable of generating the data of FIG. 3. The multiplexed assay isconfigured to perform individual assays for two targets, namely, adominant target H and a masked target L. Partition 140 contains a copyof each target, with each copy provided by a respective template 142,144, which may be single-stranded or double-stranded. (Other partitionsin the assay will contain only target H, only target L, or neithertarget.) The templates have the same sequence as one another, asindicated by regions with the same line weight, which allows bothtemplates to be bound specifically by the same pair of forward andreverse primers (F and R) and by the same probe (P). Accordingly, asignal from both target assays can be detected in the same detectionchannel. Here, the probe includes a nucleic acid 146 (anoligonucleotide) that is labeled with a luminophore 148 (e.g., afluorophore), to produce a luminescence signal. The probe also mayinclude a quencher that quenches photoluminescence from the luminophore,until the quencher and luminophore are spaced farther apart (e.g., bydegradation of the probe and/or hybridization of the probe to anamplicon, among others).

Template 144 is longer than template 142 due to the presence of anadditional sequence 150 located between the primer binding sites and notpresent in template 142. (Sequence 150 is represented by as a thickerline segment.) As a result, amplification of target H, relative totarget L, produces a shorter amplicon, which causes amplification oftarget H to proceed more efficiently than amplification of target L. Thepresence of target L is masked in the partition. Data from the assay canbe processed as described in Section II.

FIG. 7 shows a partition 160 from an exemplary multiplexed digital assaycapable of generating the data of FIG. 4. The multiplexed assay isconfigured to perform individual assays for two targets, namely, adominant or masking target L and a masked target H. Partition 160 issimilar to partition 140 of FIG. 6 except that the probe (P) has beenreplaced with a generic reporter 162, such as an intercalating dye, thatphotoluminesces more strongly when bound to double-stranded nucleic acid(e.g., amplicon). (Reporter 162 interchangeably may be termed a dye.)Light emission from dye 162 may be detected in a single optical channelor two or more optical channels.

Copies/molecules of dye 162 may bind to amplicon in a length-dependentmanner, such that more copies of the reporter bind to a copy of thelonger amplicon generated from template 144, relative to a copy of theshorter amplicon generated from template 142. As a result, amplificationof target from template 144, in the absence of template 142, generates ahigher intensity fluorescence signal than target amplification fromtemplate 142 (i.e., template 144 now corresponds to target H andtemplate 142 to target L (compare with FIG. 6)). However, the shortertarget masks the longer target, as described for FIG. 6, so target H isthe masked target in this assay configuration. (Also see Example 2 ofSection V.)

FIG. 8 shows a partition 180 from another exemplary multiplexed digitalassay capable of generating the data of FIG. 3. The partition contains acopy of a template 182 for a dominant target H and a copy of a template184 for a maskable target L. Templates 182 and 184 have the samesequence except for a sequence variation 186, such as a singlenucleotide difference, in the respective primer binding sites for thereverse primer (R). More particularly, the reverse primer forms aperfect duplex with template 182 and a mismatched duplex with template184. As a result, target H is amplified more efficiently than target L,which masks the presence of target L when both targets are present inthe same partition. The amplification inefficiency caused by themismatched reverse primer also causes a lower detected signal for targetL in the absence of target H.

FIG. 9 shows a partition 200 from yet another exemplary multiplexeddigital assay capable of generating the data of FIG. 3. The partitioncontains a copy of a template 202 for a dominant target H and a copy ofa template 204 for a maskable target L. Templates 202 and 204 have thesame sequence except for a sequence variation 206, such as a singlenucleotide polymorphism, in the respective probe binding sites for theprobe (P). More particularly, the probe forms a perfect duplex with anamplicon for target H and a mismatched duplex with an amplicon fortarget L. As a result, targets H and L are amplified with the sameefficiency, but the probe binds more efficiently to amplicon for targetH relative to the amplicon for target L, which masks the presence oftarget L, when both are present in the same partition. The bindinginefficiency caused by the mismatched probe also causes a lower signalfor target L in the absence of target H.

FIG. 10 shows a partition 220 from an exemplary multiplexed digitalassay capable of generating the data of FIG. 3. Partition 220 contains acopy of a template 222 for target H and a copy of a template 224 fortarget L. Targets H and L are composed of different sequences, asrepresented by a thicker template 222 and a thinner template 224.Accordingly, targets H and L are amplified by different pairs of primers(F_(H) and R_(H) or F_(L) and R_(L)), with amplification reported bydifferent probes (P_(H) and P_(L)). Target H is shorter than target Land is amplified substantially more efficiently, which gives target H ahigher endpoint signal and can make target H dominant to target L.

FIG. 11 shows a partition 230 from an exemplary multiplexed digitalassay performed with a single detection channel. The partition containsa copy of target A, which may be maskable, and a copy of target B, whichmay at least partially mask the presence of target A, or vice versa,according to reaction conditions (e.g., relative primer concentrationsand/or the annealing temperature for primer hybridization, amongothers). The targets are composed of different sequences, as representedby a thicker template 232 and a thinner template 234 for the targets.The targets are amplified with different respective pairs of forward andreverse primers (F_(A) and R_(A) or F_(B) and R_(B)). Amplification isreported for both targets by the same generic reporter 162 (e.g., anintercalating dye).

FIG. 11A shows a schematic graph of fluorescence intensity data that maybe collected from five sets of partitions (lanes 1-5) in a singledetection channel for the multiplexed assay of FIG. 11. Each setcontains the same concentrations of targets A and B, forward primer A,and reverse primers A and B, and a variable concentration of forwardprimer B for target B.

Lane 1 represents signal detected from generic reporter 162 binding tothe amplicon for target A only; target B exhibits no amplification inlane 1 due to the absence of forward primer B. A partition population(here, a band (A)) containing target A only has an increased signalamplitude and is well separated from the negative partitions, indicatedby ( ).

Lanes 2 and 3 represent signal detected with lower concentrations offorward primer B. The three target-positive bands (A, B, and AB) form asingle cluster, with the constituent populations not resolved from eachother. The position of the band for target A may (or may not) change(e.g., decrease in amplitude) due to competition with target B. However,the band for target A is shown here with the same position in all lanes,to simplify the presentation.

Lanes 4 and 5 show B-positive bands resolved from the band for target A.In lane 4, the two B-positive populations overlap. The data from lane 4may be utilized to determine the concentrations of both targets asdescribed above for FIG. 3. In particular, all of the lane 4 data may beutilized to calculate the concentration of target B, and only the lowesttwo bands to calculate the concentration of target A. In lane 5, allfour populations are resolved, which allows all of the target 5 data tobe used for calculation of the concentration of each target. In othercases, the concentrations of at least a pair of primers (e.g., bothF_(B) and R_(B) for target B), the annealing temperature, and/or thelike, may be adjusted to provide sufficient resolution of the A-onlypopulation from at least one other target-positive population.

FIG. 12 shows a partition 240 from an exemplary multiplexed digitalassay performed with two detection channels. The partition contains acopy of a masked target A and a copy of a masking target B. The targetsare composed of different sequences, as represented by a thickertemplate 232 and a thinner template 234 for the targets. The targets areamplified with different pairs of forward and reverse primers (F_(A) andR_(A) or F_(B) and R_(B)). Amplification is reported by different probes(P_(A) and P_(B)) labeled with distinct fluorophores 148, 242.Fluorescence from each fluorophore is detectable in a differentdetection channel as a different waveband of emitted light.

FIG. 12A shows a schematic scatter plot of amplification data collectedfrom the multiplexed assay of FIG. 12. Fluorescence intensity valuesdetected in each channel for each partition are plotted with respect toeach channel axis. Fluorophores 148 and 242 are detected substantiallyexclusively in channels 1 and 2, respectively. Fluorescence intensitydata for the partitions create three distinct clusters, which arenumbered sequentially (as 0, 1, and 2) and are identified adjacent eachcluster according to target content. Cluster “0” represents emptypartitions, namely, partitions negative for both targets, indicated as (). Cluster “1” represents partitions positive only for target A,indicated as (A). The signal intensity of partitions that are A-positiveare relatively close to those that are negative for both targets, but anaccurate partition count can be obtained from cluster 1. Cluster “2”represents partitions positive for target B, namely a population 250 ofpartitions that contain only target B, indicated as (B), and anotherpopulation 252 of partitions that contain both targets A and B,indicated as (AB). However, populations 250, 252 are not well resolvedand thus cannot provide accurate partition counts for (AB) and (B)separately. Accordingly, the concentration of target A can be calculatedusing only counts N₀ and N₁ from clusters 0 and 1, respectively, thatis, from only a portion of the data that excludes all B-positivepartitions (i.e., excluding all counts from cluster 2 (N₂):

λ_(A)=ln(N ₀ +N ₁)−ln(N ₀)  (20)

The concentration of target B, on the other hand, can be calculatedusing all of the partition counts:

λ_(B)=ln(N ₀ +N ₁ +N ₂)−ln(N ₀ +N ₁)  (21)

In other cases, each population of partitions having a particular targetcontent may be resolved from other populations with different targetcontents. For example, the AB and B populations of FIG. 12A could form apair of resolved clusters, which would allow calculation of theconcentration of each of targets A and B using data from each of thefour partition populations/clusters. More particularly, a total count ofpartitions could be based on partitions from all of thepopulations/clusters.

IV. MULTIPLEXED ASSAYS WITH A SPECIFIC REPORTER AND A GENERIC REPORTER

This section describes exemplary multiplexed assays for two or moretargets performed with at least one specific reporter and a genericreporter; see FIGS. 13-18. The assays shown here are detected in twooptical channels. However, detection may be in only one optical channelor three or more optical channels. Any of the multiplexed assaysdisclosed herein, whether performed only with one or more specificreporters, only with one or more generic reporters, or a combination ofspecific and generic reporters, may benefit from data exclusion ifpartition populations with different target content are not resolvedfrom one another in the data.

Intercalating dyes, such as SYBR® Green dye and EvaGreen® dye, arecommonly used in PCR to detect production of an amplicon. Assays basedon intercalating dyes as reporters are popular as they are lessexpensive than probe-based assays (e.g., with TaqMan® probes, molecularbeacon probes, Scorpion® primers, etc.). However, assays withintercalating dyes are not universally used for at least two reasons.First, these assays are completely reliant on the specificity of theprimers to ensure amplification of only the target of interest. For someloci, using only a pair of primers does not provide sufficientspecificity, and further specificity is achieved through a labeledprobe. The probe increases the specificity of target detection, asunwanted products that are nonspecifically amplified by the primers willnot contribute to the detected signal if the probe does not bind tothese products. Second, it can be difficult to quantify two targets of asample in the same well, when using an intercalating dye as the reporterfor both targets. However, it is generally desirable to normalize datafor a target of interest to at least one reference target (e.g., ahousekeeping gene(s)) in the sample. Without the ability to multiplex inthe same well, the sample must be split into multiple wells for separateassays to quantify the target of interest and the reference target. Thisapproach can introduce error, such as pipetting variability betweenwells. The present disclosure enables multiplexed assays with anintercalating dye as a reporter.

A generic reporter can be used, as described below, in combination withat least one specific reporter (e.g., at least one labeled probe) in adigital assay for accurate single-well quantification. This approach isattractive as it permits investigators to achieve a significant costsavings, while also reducing experimental error, because two or moretargets can be quantified simultaneously in a single well.

The multiplexed assays disclosed in this section permit researchersperforming gene expression analysis to purchase probe-based assays fortheir reference genes, and to use primer-based assays with a genericreporter for their genes of interest (GOI). If the researcher isstudying many genes, there can be a significant cost savings as only oneor a few probes need to be purchased and the other assays can bedetected with a generic reporter. Single-well quantification of a targetof interest and a reference target can be achieved, thereby eliminatingpipetting variability that currently hinders the use of an intercalatingdye for gene expression or other studies. In general, measuring morethan one target in the same multiplexed assay allows forauto-normalization of the amount of DNA loaded, as used in copy numbervariation or other applications.

FIG. 13 shows a partition 260 from another exemplary multiplexed digitalassay performed with two detection channels. The multiplexed assay isconfigured as in FIG. 11, except that the specific reporter for target Ais replaced with generic reporter 162, which is capable of binding tothe amplicons for both target A and target B. The assay has a specificreporter, probe B (P_(B)), for target B. Target A is masked, at leastpartially, by the presence of target B. However, the concentrations ofboth targets can be determined as shown below.

FIG. 14 shows a schematic scatter plot of amplification data that may becollected from the multiplexed digital assay of FIG. 13. Fluorescenceintensity values detected in each channel for each partition are plottedwith respect to each channel axis. Dye 162, which reports amplificationof both targets A and B, and fluorophore 242, which reportsamplification of target B only, are detected substantially exclusivelyin channels 1 and 2, respectively. In other words, there is sufficientspectral separation between the probe-based assay and the intercalatingdye to distinguish separate clusters. Fluorescence intensity data forthe partitions create three distinct clusters, which are numberedsequentially (as 0, 1, and 2) and are identified adjacent each clusteraccording to target content as in FIG. 12A. Cluster “2” representspartitions positive for target B, namely a population 270 of partitionsthat contain only target B, and another population 272 of partitionsthat contain both targets A and B. However, populations 270, 272 are notwell resolved and thus cannot provide accurate partition counts for (AB)and (B) separately. Accordingly, the concentrations of targets A and Bcan be calculated using Equations 20 and 21, respectively. In anexemplary embodiment, target A is RPP30 and target B is MRG.

In other cases, each population of partitions having a particular targetcontent may be resolved from other populations with different targetcontents. For example, the AB and B populations of FIG. 14 (or FIG. 14Abelow) could form a pair of resolved clusters, which would allowcalculation of the concentration of each of targets A and B using datafrom each of the four partition populations/clusters. More particularly,a same total count of partitions used for calculating each concentrationcould be based on partitions from all of the populations/clusters.

FIG. 14A shows a schematic view of a partition 274 from an exemplarymultiplexed digital assay performed with a single detection channel. Themultiplexed assay is configured as in FIG. 13, except that probe B(P_(B)) for target B is labeled with fluorophore 148, which allowsfluorescence from generic reporter 162 and fluorophore 148 of probe B tobe detected in the same channel. Partition 274 contains a copy of targetA and target B, each amplifiable with a different pair of forward andreverse primers (F_(A) and R_(A) or F_(B) and R_(B)).

FIG. 14B shows a schematic graph of fluorescence intensity datacollected from five sets of partitions (lanes 1-5) in a single detectionchannel for the multiplexed digital assay of FIG. 14A. Partition sets1-5 contain the same concentrations of targets A and B and of both pairsof primers, but have a variable concentration of the probe (P_(B))specific for target B.

Lane 1 represents signal detected from generic reporter 162 only (probeB is not present). Partition populations containing target A only (A),target B only (B), or both targets A and B (AB) have an increased signalamplitude and are well separated from the negative partitions ( ).However, the three target-positive populations are not resolved fromeach other.

Lanes 2-5 represent signal detected from the assay performed with anincreasing concentration of probe B. The position of the band for targetA may not change, but each of the B-positive bands migrates upward inthe graph with increased concentrations of probe B, to produce a clearseparation of the B-negative and B-positive populations. Accordingly,the data of lane 4 or lane 5 can be utilized to accurately determine theconcentrations of targets A and B. The concentration of target A can bedetermined from only a subset of the data, namely, from partition countsof the negative band ( ) and the A-only band (A), with Equation 20. Theconcentration of target B can be determined from partition counts of allof the bands in either lane, with Equation 21.

FIG. 15 is a schematic view of a partition 280 from an exemplarymultiplexed digital assay performed with two detection channels andthree unlinked targets (A-C). The multiplexed assay configuration ofFIG. 15 is related to that of FIG. 13, except that an assay for another,unrelated target (C) is performed (a copy of a template 282 for thetarget is present in partition 280). Target C is amplified with targetC-specific primers (F_(C) and Rc), and amplification is reported by atarget C-specific probe (P_(C)), which is labeled with fluorophore 148.The exact number and arrangement of clusters will depend on variousparameters, such as the condition(s) limiting each amplificationreaction, the channels used for detection, the effective concentrationof each reporter, the spectral characteristics of each reporter, and thelike. Any suitable combination of the parameters can be adjusted toachieve a cluster arrangement that permits accurate concentrationdetermination.

FIG. 16 is a schematic scatter plot of amplification data collected fromthe multiplexed digital assay of FIG. 15. Fluorescence intensity (FL.INT.) values detected in each channel for each partition are plottedwith respect to each channel axis. Dye 162 and fluorophore 148 (targetC) are detected substantially exclusively in channel 1, and fluorophore242 (target B) is detected substantially exclusively in channel 2.Fluorescence intensity data for the partitions create six distinctclusters, which are numbered sequentially (as 0, 1, 2, 3, 4, and 5) andare identified adjacent each cluster in parentheses according to targetcontent. Cluster 2 represents partitions positive for target C, namely apopulation 290 of partitions that contain only target C, and anotherpopulation 292 of partitions that contain both targets A and C. Cluster5 represents partitions positive for targets B and C, namely, apopulation 294 of partitions that contain only target B and C, andanother population 296 of partitions that contain all three targets.However, populations 290, 292 and populations 294, 296 are not wellresolved and thus cannot provide accurate partition counts for eachindividual population. However the concentrations of targets A, B, and Ccan be calculated with proper combination and exclusion of clustercounts.

The concentration of target A can be calculated as follows:

λ_(A)=ln(N ₀ +N ₁)−ln(N ₀)  (22)

Alternatively, the concentration of target A can be calculated byincluding clusters 3 and 4 and excluding all C-positive clusters:

λ_(A)=ln(N ₀ +N ₁ +N ₃ +N ₄)−ln(N ₀ +N ₃)  (23)

The concentration of target B can be calculated using counts from allthe partitions as follows:

λ_(B)=ln(N _(tot))−ln(N ₀ +N ₁ +N ₂)  (24)

Alternatively, the concentration of target B can be calculated byexcluding counts from all C-positive clusters (i.e., by excluding N₂ andN₅):

λ_(B)=ln(N ₀ +N ₁ +N ₃ +N ₄)−ln(N ₀ +N ₁)  (25)

Clusters 2 and 5 can be excluded from the calculation because each isexpected to have the same fraction of A-positive partitions as all ofthe partitions considered collectively. However, generally, it isdesirable to use more of the data for a calculation, if accuracy doesnot suffer.

The concentration of target C can be calculated by excluding allB-positive clusters (i.e., by excluding N₃, N₄, and N₅):

λ_(C)=ln(N ₀ +N ₁ +N ₂)−ln(N ₀ +N ₁)  (26)

Alternatively, the concentration of target C can be calculated usingcounts from all the partitions as follows:

λ_(C)=ln(N _(tot))−ln(N ₀ +N ₁ +N ₃ +N ₄)  (27)

FIG. 17 shows another schematic scatter plot of amplification datacollected from the multiplexed digital assay of FIG. 15 and plotted asin FIG. 16, except that cluster 3 (B) and cluster 4 (AB) have merged toform a new cluster 3, and cluster 5 has been renamed as cluster 4.

Here, again, the concentrations of all targets can be calculated. Theconcentration of target A can be calculated with Equation 22. Theconcentration of target B can be calculated using counts from all of thepartitions with Equation 24. Alternatively, the concentration of targetB can be calculated by excluding counts from all C-positive clusters(i.e., by excluding N₂ and N₄):

λ_(B)=ln(N ₀ +N ₁ +N ₃)−ln(N ₀ +N ₁)  (28)

The concentration of target C can be calculated by excluding allB-positive clusters (i.e., by excluding N₃ and N₄) using Equation 26.Alternatively, the concentration of target C can be calculated usingcounts from all the partitions as follows:

λ_(C)=ln(N _(tot))−ln(N ₀ +N ₁ +N ₃)  (29)

FIG. 18 shows yet another schematic scatter plot of amplification datacollected from the multiplexed digital assay of FIG. 15. Fluorescenceintensity (FL. INT.) values detected in each channel for each partitionare plotted with respect to each channel axis. Generic reporter 162 isdetected in both channels 1 and 2, fluorophore 148 (target C)substantially exclusively in channel 1, and fluorophore 242 (target B)substantially exclusively in channel 2. Fluorescence intensity data forthe partitions create five distinct clusters, which are numberedsequentially (as 0, 1, 2, 3, and 4) and are identified adjacent eachcluster in parentheses according to target content. The clusters of FIG.18 correspond to those of FIG. 17; the strategies of FIG. 17 forcalculating concentration of targets A-C can be applied to the clustercounts of FIG. 18.

V. EXAMPLES

This section describes selected aspects and embodiments of the presentdisclosure related to multiplexed digital assays. The assays may involveany combination of target masking, data exclusion for calculation of atleast one target level, use of the same generic reporter for at leasttwo targets in the same multiplexed assay, and/or use of a specificreporter and a generic reporter for the same multiplexed assay, amongothers. These examples are intended for illustration only and should notlimit or define the entire scope of the present disclosure.

Example 1 Analysis of Library Quality

This example describes use of a multiplexed digital assay in droplets tomeasure levels of empty and filled members of a library, where thepresence of an empty library member in a droplet masks the presence of afilled member in the same droplet; see FIGS. 19 and 20.

FIG. 19 shows a schematic view of a partition 300 from an exemplarymultiplexed digital assay to determine the quality of a library, such asan Ion Torrent® library. The library may be a Next-Generation Sequencing(NGS) library constructed in vitro by attaching a copy of an adapter(“A”) to each end of a population of fragments of interest. Members ofthe library can include an empty member 302 produced by ligation of apair of adapter copies to each other in an inverted orientation, withoutan intervening insert. Empty members 302 do not provide useful sequenceinformation and thus are an undesired component of the library. The goalof library construction is to generate a high-complexity population offilled members, such as member 304, each containing one of the fragmentsof interest as an insert flanked by attached copies of the adapter.

The quality of an NGS library can be characterized in a multiplexeddigital assay. Members of the library can be distributed at partialoccupancy to partitions, such as droplets. The droplets can contain oneor more primers (F) to amplify members of the library and anadapter-specific probe (P) to report amplification. Here, the sameprimer (F) can function as a forward primer and a reverse primer foramplification by binding specifically to both copies of the adapter.

FIG. 20 shows a graph of amplification data collected from themultiplexed assay of FIG. 19 performed in droplets. The probe is labeledwith FAM dye. Fluorescence intensity measured from the FAM dye isplotted as a function of event (or droplet) number. The empty librarymember (A-A), which is very short, amplifies very efficiently andproduces a stronger signal than a filled library member, which masks thesignal from the filled library member, if present in the same droplet.Accordingly, the concentrations of empty and filled library members canbe calculated as described for FIG. 3. In particular, the concentrationof empty library members, which are dominant/masking, can be calculatedusing counts from all the droplets, and the concentration of the filledlibrary members, which are masked, can be calculated by excluding thecounts for the empty library members. Exemplary concentrations obtainedwith this approach are shown below for four libraries:

Library Total Concentration A-A Concentration N13 574 19.7 N14 539 30.0N15 631 29.8 N16 581 9.17Further aspects of library construction and characterization by digitalassay are described in U.S. patent application Ser. No. 13/562,198,filed Jul. 30, 2012, which is incorporated herein by reference.

Example 2 Measurement of Spliced and Unspliced RNA Levels

This example describes use of a multiplexed digital assay in droplets tomeasure levels of unspliced and spliced transcripts in a sample; seeFIGS. 21-24.

FIG. 21 shows a schematic view of a partition 320 from an exemplarymultiplexed digital assay to quantify spliced and unspliced species in asample. Partition 320 contains a copy of an unspliced species(Exon-Intron-Exon) and a spliced species (Exon-Exon). Each species isamplifiable with the same pair of forward and reverse primers (F and R).Amplification is reported via light emission from dye 162, which bindsnonspecifically to amplicons generated by amplification of target fromboth templates.

FIG. 22 shows a graph of amplification data collected from themultiplexed assay of FIG. 21 performed in droplets with primers to exons7 and 8 of GAPDH, and SBYR® Green as dye 162. The expected ampliconsized is 64 bp for the spliced species with the intron removed.Additional droplets with higher intensity signals are visible with anamplicon size of 257 bp, which may be produced from immature mRNAmolecules where the introns have not yet been excised. The 64 bp speciesoutcompetes the 257 bp species and masks its signal. Note the inversionof signal since the dye is intercalating and produces signalproportional to amplicon length (also see FIGS. 4 and 7). Theconcentration calculated for the 257 bp species was confirmed in anassay using intron primers to amplify the intron, without competitionfrom the shorter, spliced species.

FIG. 23 shows a partition 340 from another exemplary multiplexed digitalassay to quantify spliced and unspliced species in a sample. Thepartition contains a copy of an unspliced species 342(Exon-Intron-X-Intron-Exon) and two distinct spliced species 344, 346(Exon-X-Exon and Exon-Exon). Each species is amplifiable with the samepair of forward and reverse primers (F and R). Amplification of eachspecies is reported by the same dye 162.

FIG. 24 shows a graph of amplification data collected from themultiplexed assay of FIG. 23 performed in droplets with primers to exons3 and 4 of actin. Here, three amplicons of different length (76 bp, 117bp, and 517 bp) are visible. The shortest amplicon outcompetes bothlonger amplicons, and the intermediate amplicon outcompetes the longestamplicon. Concentrations of the three targets can be computed asdescribed in Section II.

Further aspects of multiplexed assays for insertions and deletions aredescribed in U.S. Provisional Patent Application Ser. No. 61/692,635,filed Aug. 23, 2012, which is incorporated herein by reference.

Example 3 Multiplexed Digital Assay with a Reporter-less Target

This example describes a multiplexed digital assay performed inpartitions to measure levels of a pair of targets, and using a singlereporter that binds to only one of the targets; see FIGS. 25 and 26.

FIG. 25 shows an exemplary partition 360 from an exemplary multiplexeddigital assay for a first target (target A) and a second target (targetB). Partition 360 may be a member of a set of partitions having copiesof target A and target B randomly distributed among the partitions. Onlya subset of the partitions each contain at least one copy of target A.Only a different subset of the partitions each contain at least one copyof target B. Yet another different subset of the partitions each containat least one copy of target A and at least one copy of target B, withpartition 360 shown in FIG. 25 as an example.

Each partition may include or be an isolated portion of a mixture. Themixture may contain targets A and B (e.g., provided by templates 362 and364, respectively), and also may contain reagents sufficient to amplifyeach target under appropriate environmental conditions or manipulation(e.g., heating, cyclical heating and cooling, or the like).

Any suitable reagents may be included in the mixture and thus in eachpartition 360. The reagents may include forward and reverse primers(F_(A) and R_(A)) for amplification of target A, and forward and reverseprimers (F_(B) and R_(B)) for amplification of target B. The reagentsalso may include a specific reporter (e.g., a probe, P_(A), including alabel, photoluminophore 366). The reporter binds specifically to targetA (and/or to a region of an amplicon formed by amplification of targetA). However, the mixture does not contain a reporter that binds totarget B or that is sensitive to amplification of target B, exceptindirectly. More particularly, the presence of target B is detectablevia probe P_(A), as a change in the probe signal, only when both targetsA and B are amplified in the same partition, such as in partition 360that contains a copy of each target. Without the presence of target A ina partition, target B remains invisible and undetectable for thepartition.

FIG. 26 shows a schematic scatter plot of amplification data(fluorescence intensity) that may be collected from partitions in themultiplexed digital assay of FIG. 25. The plot shows data collected intwo channels, but photoluminophore 366 of probe A is detectable only inchannel 1. In other cases, the data may be plotted as intensity withrespect to time or event number. In any event, three clusters ofpartitions are distinguishable and are numbered sequentially (0, 1, and2) and are identified adjacent each cluster according to target content.

Cluster 0 has the lowest intensity in channel 1, because each partitionof cluster 0 is negative for target A. Since target A is absent,amplification of unprobed target B goes undetected. Therefore, target Bis distributed among the partitions of cluster 0 at about the sameaverage concentration (copy per droplet) as for the partitions in all ofthe clusters, considered collectively. Any given partition of cluster 0can be positive or negative for target B.

Cluster 2 has the highest intensity in channel 1 and contains target Aonly. Amplification of target A occurred in the partitions of cluster 2without competition from target B, which results in more efficientamplification of target A and thus a stronger signal from the probe.

Cluster 1 has an intermediate intensity in channel 1. Each partition ofcluster 1 amplifies target A and target B. The A and B amplificationreactions compete with one another for limiting reagents, which resultsin less efficient amplification of target A relative to the partitionsof cluster 2. As a result, the probe produces a signal of lesserintensity for cluster 1 than for cluster 2.

N₀, N₁, and N₂ are the number of partitions (the partition counts) inclusters 0, 1, and 2, respectively. The concentration (copies perdroplets) of target A then can be calculated as follows:

λ_(A)=ln(N ₀ +N ₁ +N ₂)−ln(N ₀)  (30)

In Equation 30, the total number of partitions is provided by summingthe partition numbers for all three clusters, and the number ofpartitions negative for target A only by the number of partitions incluster 0 (also see Equation 5 above).

The concentration (copies per droplets) of target B can be calculated asfollows:

λ_(B)=ln(N ₁ +N ₂)−ln(N ₂)  (31)

In Equation 31, the total number of partitions may be determined bysumming the partition numbers for clusters 1 and 2, and the number ofpartitions negative for target B by the number of partitions in cluster2 (also see Equation 5 above). Cluster 0 can be excluded from thecalculation of the target B concentration because the concentration oftarget B in cluster 0 is expected to be the same (absent statisticalerror) as in clusters 1 and 2, taken collectively.

Example 4 Selected Embodiments I

This example describes selected embodiments of multiplexed digitalassays, presented as a series of numbered paragraphs.

1. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing partitions containing a first target and asecond target at partial occupancy, with a plurality of the partitionseach containing a copy of both targets; (B) performing a reactioncorresponding to each target in the partitions; (C) collecting data fromat least one reporter for each reaction; and (D) determining aconcentration of the first target from only a subset of the data thatselectively excludes partitions positive for the second target.

2. The method of paragraph 1, wherein the partitions are droplets,wherein the step of collecting data is performed by a detectionassembly, wherein the step of determining is performed with a dataprocessor, or any combination thereof.

3. The method of paragraph 2, wherein the step of providing partitionsincludes a step of forming droplets containing the first target and thesecond target.

4. The method of any of paragraphs 1 to 3, wherein the step of providingpartitions includes a step of forming partitions including a samplecontaining the first target and the second target.

5. The method of any of paragraphs 1 to 4, wherein the first target andthe second target each include a sequence of nucleotides.

6. The method of any of paragraphs 1 to 5, wherein a plurality of thepartitions each contain a single copy of the first target and aplurality of the partitions each contain a single copy of the secondtarget when the partitions are provided.

7. The method of any of paragraphs 1 to 6, wherein the reactioncorresponding to each target amplifies nucleic acid.

8. The method of any of paragraphs 1 to 7, wherein the reactioncorresponding to each target is catalyzed by at least one enzyme.

9. The method of paragraph 8, wherein the at least one enzyme includes apolymerase, a ligase, or both.

10. The method of paragraph 8, wherein the at least one enzyme includesan enzyme connected to an oligonucleotide that binds specifically to oneof the targets.

11. The method of any of paragraphs 1 to 10, wherein the at least onereporter for each reaction includes a reporter that binds to nucleicacid.

12. The method of paragraph 11, wherein a reporter for at least one ofthe reactions is a labeled probe that includes an oligonucleotide.

13. The method of paragraph 12, wherein the step of performing areaction produces a first amplicon corresponding to the first target anda second amplicon corresponding to the second target, wherein a reporterfor the reaction corresponding to the first target binds to the firstamplicon and the second amplicon.

14. The method of paragraph 13, wherein a reporter for the reactioncorresponding to the second target is a specific reporter that bindsspecifically to the second amplicon relative to the first amplicon, orwherein the data is collected from only one reporter.

15. A method of performing a multiplexed digital assay, the methodcomprising: (A) amplifying a first target and a second target inpartitions containing each target at partial occupancy; (B) collectingdata for amplification of the first target and the second target, thesecond target at least partially masking the presence of the firsttarget in the data for a plurality of the partitions that contain a copyof each target; and (C) determining a concentration of the first targetfrom only a subset of the data that selectively excludes partitionspositive for the second target.

16. The method of paragraph 15, wherein the partitions are droplets.

17. The method of paragraph 15 or 16, wherein each of the targets ispresent in the partitions at an average of less than five copies perpartition.

18. The method of any of paragraphs 15 to 17, wherein each of aplurality of the partitions contains no copy of either target.

19. The method of any of paragraphs 15 to 18, wherein the concentrationis an average number of copies per partition of the first target.

20. The method of any of paragraphs 15 to 19, wherein the step ofdetermining is based on a value for a total number of partitionsrepresented by the subset of the data.

21. The method of any of paragraphs 15 to 20, wherein the step ofdetermining is based on a value representing a total number ofpartitions in the subset of the data that are negative for the firsttarget, or on a value representing a total number of partitions in thesubset of the data that are positive for the first target.

22. The method of any of paragraphs 15 to 21, wherein the step ofdetermining is based on Poisson statistics.

23. The method of any of paragraphs 15 to 22, further comprising a stepof determining a concentration of the second target based on at least asubset of the data.

24. The method of paragraph 23, wherein the steps of amplifying andcollecting are performed for a third target, and wherein theconcentration of the second target is determined from only a subset ofthe data that excludes partitions positive for the third target.

25. The method of any of paragraphs 15 to 24, wherein amplification ofthe second target suppresses amplification of the first target inpartitions containing a copy of each target.

26. The method of any of paragraphs 15 to 25, wherein at least one sameprimer is used for amplification of the first target and the secondtarget.

27. The method of paragraph 26, wherein a same pair of forward andreverse primers is used for amplification of the first target and thesecond target.

28. The method of paragraph 26 or 27, wherein amplification of the firsttarget and the second target generates respective first and secondamplicons, and wherein the first amplicon is substantially longer thanthe second amplicon, or vice versa.

29. The method of paragraph 28, wherein the second amplicon is at least25%, 50%, or 100% longer than the first amplicon, or vice versa.

30. The method of any of paragraphs 15 to 29, further comprising a stepof identifying from the data at least one cluster of partitions positivefor the second target, wherein the step of determining is performed withthe at least one cluster excluded.

31. The method of paragraph 30, further comprising a step of generatingat least one plot of the data, wherein the step of identifying isperformed based on the least one plot.

32. The method of any of paragraphs 15 to 31, further comprising a stepof identifying the subset of the data based on comparing one or moresignal values for each partition with one or threshold values or ranges.

33. The method of paragraph 32, wherein the one or more signal valuesrepresent a detected intensity of emitted light.

34. The method of any of paragraphs 15 to 33, wherein amplification ofthe first target and second target generates respective first and secondamplicons, and wherein the each of the partitions contains a samereporter that binds to the first amplicon and the second amplicon.

35. The method of paragraph 34, wherein the same reporter is a dye thatbinds to double-stranded nucleic acid.

36. The method of paragraph 34 or 35, wherein the same reporter binds tothe first and second amplicons without substantial sequence specificity.

37. The method of paragraph 34, wherein the same reporter is a probethat binds specifically to the first and second amplicons.

38. The method of any of paragraphs 15 to 37, wherein the step ofamplifying also amplifies a third target, wherein the data collected isalso for amplification of the third target in the partitions, whereinthe third target at least partially masks the presence of the secondtarget for partitions that contain a copy of each of the second andthird targets, further comprising a step of determining a concentrationof the second target from only a subset of the data that selectivelyexcludes partitions positive for the third target.

39. The method of paragraph 38, wherein the third target at leastpartially masks the presence of the first target in the data forpartitions that contain a copy of each of the first and third targets.

40. The method of paragraph 38 or 39, wherein the step of determining aconcentration of the first target is performed with only a subset of thedata that selectively excludes partitions positive for the second targetand selectively excludes partitions positive for the third target.

41. The method of any of paragraphs 15 to 40, wherein the subset of thedata from which the concentration of the first target is determinedexcludes at least a majority of the partitions positive for the secondtarget.

42. The method of paragraph 41, wherein the subset of the data excludesall partitions that test positive for the second target.

43. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing partitions containing a first target and asecond target at partial occupancy, with a plurality of the partitionseach containing a copy of both targets; (B) amplifying the targets inthe partitions; (C) collecting data for amplification of the firsttarget and the second target, the data being plottable to produce threeor more clusters of data points with each cluster representing at leastone different combination of the targets in the partitions; and (D)determining a concentration of the first target from only a subset ofthe data that excludes one or more of the clusters.

44. The method of paragraph 43, wherein one of the three or moreclusters represents partitions that are negative for both targets.

45. The method of paragraph 43 or 44, wherein one of the three or moreclusters represents partitions that are positive for the first targetand negative for the second target.

46. The method of any of paragraphs 43 to 45, wherein one of the threeor more clusters represents partitions that are positive for the secondtarget without regard to whether the partitions are positive or negativefor the first target.

47. The method of any of paragraphs 43 to 46, further comprising a stepof plotting the data to form the clusters.

48. The method of paragraph 47, wherein the data is plotted asluminescence intensity as a function of event number.

49. The method of paragraph 47, wherein the data is plotted to generatea two-dimensional scatter plot having a first axis and a second axis,and wherein each axis represents intensity for a different waveband oflight.

50. The method of any of paragraphs 43 to 49, wherein one of the threeor more clusters includes a pair of sub-clusters representing differentcombinations of the targets, and wherein the sub-clusters overlap oneanother.

Example 5 Selected Embodiments II

This example describes selected embodiments of multiplexed digitalassays, presented as a series of numbered paragraphs.

1. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing a mixture including a plurality of targets andreagents sufficient for amplification of each of the targets, whereinthe plurality of targets includes a first target; (B) forming partitionseach including a portion of the mixture; (C) amplifying the plurality oftargets in the partitions; (D) collecting data from a plurality of thepartitions for amplification of each of the plurality of targets; (E)identifying partition populations in the data, each partition populationhaving a different target content of the plurality of targets; and (F)calculating a level of the first target from only a portion of the data,wherein the portion of the data is collected from only a subset of theplurality of the partitions, wherein the subset of partitions excludestwo or more of the partition populations, and wherein the two or moreexcluded partition populations include at least a pair of excludedpartition populations that are not resolved from each other in the data.

2. The method of paragraph 1, wherein the at least a pair of excludedpartition populations that are not resolved from each other include afirst partition population that is positive for the first target and asecond partition population that is negative for the first target.

3. The method of paragraph 1 or 2, wherein the step of identifyingpartition populations includes a step of plotting the data to form threeor more clusters of data points, wherein each cluster represents atleast one of the partition populations, wherein the portion of the dataexcludes one or more of the clusters, and wherein the one or moreexcluded clusters include a heterogeneous cluster composed of at leasttwo overlapping sets of data points representing partition populationshaving different target content from one another.

4. The method of paragraph 3, wherein each data point represents one ormore intensity values of light detected from each partition of theplurality of the partitions, and wherein the step of plotting the dataincludes a step of plotting the data with respect to at least oneintensity axis.

5. The method of paragraph 4, where the step of plotting the dataincludes a step of plotting the data with respect to at least twointensity axes.

6. The method of any of paragraphs 1 to 5, wherein the step ofcollecting data includes a step of obtaining at least one value for eachpartition of the plurality of the partitions, and wherein the step ofidentifying is based on comparing the at least one signal value to atleast one threshold or range for each partition of the plurality of thepartitions.

7. The method of any of paragraphs 1 to 6, wherein the plurality oftargets includes a second target, and wherein the portion of the dataexcludes each partition testing positive for the second target.

8. The method of paragraph 7, wherein the step of amplifying isperformed with a same pair of primers for the first target and thesecond target.

9. The method of paragraph 7 or 8, wherein a subset of the plurality ofthe partitions each contain at least one copy of the first target and atleast one copy of the second target, and wherein amplification of thesecond target reduces amplification of the first target in partitionseach containing a copy of the first target and the second target.

10. The method of any of paragraphs 7 to 9, wherein the plurality oftargets includes a third target, further comprising a step ofcalculating a level of the second target based on only a portion of thedata that excludes partitions testing positive for the third target.

11. The method of any of paragraphs 1 to 10, wherein the level is anaverage number of copies per partition of the first target in the subsetof partitions, and wherein the average number is expected to be thesame, absent statistical error, as an average number of copies perpartition of the first target in the two or more excluded partitionpopulations taken collectively.

12. The method of any of paragraphs 1 to 11, wherein the step ofidentifying partition populations includes a step of identifying apartition population that tests negative for each of the plurality oftargets.

13. The method of paragraph 12, wherein the subset of partitions doesnot include the partition population that tests negative for each of theplurality of targets.

14. The method of any of paragraphs 1 to 13, wherein the step ofidentifying partition populations includes a step of identifying atleast a pair of the partition populations that are not resolved fromeach other in a plot of the data.

15. The method of paragraph 14, wherein the subset of partitions doesnot include the at least a pair of partition populations that are notresolved from each other in a plot of the data.

16. The method of any of paragraphs 1 to 15, wherein only a first subsetof the plurality of the partitions each contain at least one copy of thefirst target, wherein only a distinct second subset of the plurality ofthe partitions each contain at least one copy of a second target of theplurality of targets, and wherein a third subset of the plurality of thepartitions each contain at least one copy of the first target and thesecond target.

17. The method of any of paragraphs 1 to 16, further comprising a stepof obtaining a first partition count and a second partition count fromthe subset of partitions, wherein the step of calculating a level isperformed with the first partition count and the second partition count.

18. The method of paragraph 17, wherein the first partition count is atotal number of partitions in the subset of partitions, and wherein thesecond partition count is a number of partitions positive for the firsttarget in the subset of partitions or a number of partitions negativefor the first target in the subset of partitions.

19. The method of any of paragraphs 1 to 18, wherein the step of formingpartitions includes a step of forming droplets.

20. The method of any of paragraphs 1 to 19, wherein each target of theplurality of targets includes a sequence of nucleotides.

21. The method of any of paragraphs 1 to 20, wherein the mixtureincludes at least one reporter, and wherein the step of collecting dataincludes a step of detecting light from the at least one reporter.

22. The method of paragraph 21, wherein the at least one reporterincludes a same reporter that is sensitive to amplification of the firsttarget and a second target of the plurality of targets.

23. The method of paragraph 22, wherein the same reporter is a genericreporter.

24. The method of paragraph 23, wherein the generic reporter binds todouble-stranded nucleic acid without substantial sequence specificity.

25. The method of any of paragraphs 21 to 24, wherein the at least onereporter includes a labeled probe that bind specifically to only thefirst target of the plurality of targets.

26. The method of any of paragraphs 1 to 25, wherein the step ofcalculating a level includes a step of calculating a concentration ofthe first target.

27. The method of paragraph 26, wherein the step of calculating aconcentration of the first target includes a step of calculating anaverage number of copies of the first target per partition.

28. The method of any of paragraphs 1 to 27, wherein the step ofcalculating a level is based on Poisson statistics.

29. The method of any of paragraphs 1 to 28, wherein the step ofamplifying the first target generates an amplicon, and wherein thepartitions do not include a reporter that binds to the first target, theamplicon, or both the first target and the amplicon.

30. The method of any of paragraphs 1 to 29, wherein the plurality oftargets include a second target, and wherein the partitions include areporter that is differentially sensitive to amplification of the secondtarget in a partition relative to amplification of the first target andthe second target in the partition.

31. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing a mixture including a plurality of targets andreagents sufficient for amplification of each of the targets, whereinthe plurality of targets includes a first target and a second target;(B) forming partitions each including a portion of the mixture; (C)amplifying the plurality of targets in the partitions; (D) collectingdata from a plurality of the partitions for amplification of each of theplurality of targets, the data being plottable to form three or moreclusters of data points, with each data point representing a partition,and with each cluster having a different content of the plurality oftargets; and (E) calculating a level of the first target from only aportion of the data that excludes at least one of the clusters of datapoints, wherein the at least one excluded cluster includes aheterogeneous excluded cluster composed of overlapping sets of datapoints representing at least two excluded partition populations, whereinat least one of the excluded partition populations is positive for thefirst target, and wherein at least one of the excluded partitionpopulations is negative for the first target.

32. The method of paragraph 31, further comprising a step of plottingthe data to form the three or more clusters of data points.

33. The method of paragraph 31 or 32, further comprising a step ofidentifying each cluster according to the target content of the cluster.

34. The method of paragraph 33, wherein the portion of the data excludesat least a pair of heterogeneous clusters.

35. The method of any of paragraphs 31 to 34, wherein the plurality oftargets includes a second target, and wherein the portion of the dataexcludes data points for partitions testing positive for the secondtarget.

36. The method of any of paragraphs 31 to 35, wherein the portion of thedata excludes data points for partitions testing negative for each ofthe plurality of targets.

37. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing a mixture including a first target, a secondtarget, and reagents sufficient for amplification of each of thetargets; (B) forming partitions each including a portion of the mixture;(C) amplifying the first target and the second target in the partitions;(D) collecting data from a plurality of the partitions for amplificationof the first target and the second target; and (E) calculating a levelof the first target from only a portion of the data that excludes eachpartition testing positive for the second target.

38. The method of paragraph 37, wherein first target is not linked tothe second target in the mixture and is not linked to the second targetwhen the step of forming partitions is performed.

39. The method of paragraph 37 or 38, wherein the data is collected fromat least one reporter including a labeled probe.

40. The method of paragraph 39, wherein the labeled probe has an intactform and one or more degraded forms, and wherein the step of collectingdata includes a step of collecting at least part of the data from one ormore degraded forms of the probe.

41. The method of any of paragraphs 37 to 40, wherein a presence of thesecond target at least partially masks a presence of the first target inthe data for a plurality of the partitions each containing at least onecopy of each target.

42. The method of any of paragraphs 37 to 41, wherein the portion of thedata is collected from only a subset of the plurality of the partitions,wherein the step of calculating a level is based on a total numberpartitions in the subset of the partitions and on a number of partitionsin the subset of the partitions that are negative for the first targetor that are positive for the first target.

43. A method of performing a multiplexed digital assay, the methodcomprising: (A) amplifying a first target and a second target inpartitions containing each target at partial occupancy; (B) collectingdata for amplification of the first target and the second target, thesecond target at least partially masking a presence of the first targetin the data for a plurality of the partitions that each contain at leastone copy of each target; and (C) calculating a concentration of thefirst target from only a portion of the data that selectively excludespartitions testing positive for the second target.

44. The method of paragraph 43, wherein the step of calculating is basedon a value for a total number of partitions represented by the portionof the data.

45. The method of paragraph 44, wherein the step of calculating is basedon a value for a number of partitions in the portion of the data thattest negative for the first target, or on a value for a number ofpartitions in the portion of the data that test positive for the firsttarget.

46. The method of any of paragraphs 43 to 45, further comprising a stepof calculating a level of the second target based on at least a portionof the data.

47. The method of any of paragraphs 43 to 46, wherein the steps ofamplifying and collecting data are performed for a third target, furthercomprising a step of calculating a level of the second target from onlya portion of the data that excludes partitions testing positive for thethird target.

48. The method of any of paragraphs 43 to 47, wherein amplification ofthe second target reduces amplification of the first target inpartitions each containing a copy of each target.

49. The method of any of paragraphs 43 to 48, further comprising a stepof identifying from the data at least one cluster of partitions positivefor the second target, wherein the step of calculating is performed withthe at least one cluster excluded.

50. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing a mixture including a plurality of targets andreagents sufficient for amplification of each of the targets, whereinthe plurality of targets includes a first target and a second target;(B) forming partitions each including a portion of the mixture, withonly a first subset of the partitions containing the first target andonly a distinct second subset of the partitions containing the secondtarget, and with a distinct third subset of the partitions eachcontaining at least one copy of each of the first target and the secondtarget; (C) amplifying the plurality of targets in the partitions; (D)collecting data from a plurality of the partitions for amplification ofeach of the plurality of targets, the data being plottable to producethree or more clusters of data points with each cluster representing atleast one target combination not represented by any other of the threeor more clusters; and (E) calculating a level of the first target fromonly a portion of the data that excludes one or more of the clusters,wherein the one or more excluded clusters include an excluded clusterthat is heterogeneous with respect to the first target (i.e., only asubset the cluster is positive for the first target).

51. The method of paragraph 50, wherein the at least one targetcombination includes a combination in which each of the plurality oftargets is absent.

52. The method of paragraph 50 or 51, wherein the at least one targetcombination includes a combination in which only one target of theplurality of targets is present.

53. The method of any of paragraphs 50 to 52, wherein the at least onetarget combination includes a combination in which the first target isabsent and the second target is present.

54. The method of any of paragraphs 50 to 53, wherein every clustercontaining the second target is excluded from the step of calculating alevel.

55. The method of any of paragraphs 50 to 54, wherein the level is aconcentration, wherein the concentration is an average number of copiesper partition, and wherein the average number is expected to be thesame, absent statistical error, as an average number of copies perpartition of the first target in the one or more excluded clusters takencollectively.

56. The method of any of paragraphs 50 to 55, wherein the step ofcalculating a level uses a value for a total number of partitionsrepresented by the portion of data, and wherein partitions of the one ormore clusters are not included in the total number of partitions.

57. The method of any of paragraphs 50 to 56, wherein one of the threeor more clusters includes data points for partitions that are negativefor the first target and the second target.

58. The method of any of paragraphs 50 to 57, wherein one of the threeor more clusters includes data points for partitions that are positivefor the first target and negative for the second target.

59. The method of any of paragraphs 50 to 58, wherein the excludedcluster that is heterogeneous with respect to the first target includesdata points for partitions that are positive for the first target anddata points for partitions that are negative for the first target.

60. The method of any of paragraphs 50 to 59, further comprising a stepof plotting the data to form the clusters.

61. The method of paragraph 60, wherein the data is plotted asphotoluminescence intensity as a function of event number.

62. The method of paragraph 60, wherein the data is plotted to generatea two-dimensional scatter plot having a first axis and a second axis,and wherein each axis represents detected intensity for a differentwavelength or waveband of light.

Example 6 Selected Embodiments III

This example describes selected embodiments of multiplexed digitalassays performed with a generic reporter, presented as a series ofnumbered paragraphs.

1. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing partitions containing a first target and asecond target at partial occupancy; (B) amplifying the targets in thepartitions; (C) collecting data for amplification of the first targetand the second target in the partitions from one or more reportersincluding a generic reporter that binds to amplicons representing bothtargets; and (D) determining a concentration of the first target and thesecond target based on the data.

2. The method of paragraph 1, wherein a plurality of the partitions eachcontain a copy of both targets.

3. The method of paragraph 1 or 2, wherein the concentration of thefirst target is determined from only a subset of the data thatselectively excludes partitions positive for the second target.

4. The method of any of paragraphs 1 to 3, wherein at least one sameprimer is used for amplification of both targets.

5. The method of paragraph 4, wherein amplification of each target isperformed with a same pair of primers.

6. The method of any of paragraphs 1 to 3, wherein amplification of eachtarget is performed with a pair of primers, only one of which or neitherof which binds to the other target.

7. The method of any of paragraphs 1 to 6, wherein the generic reporterincludes a dye that binds to double-stranded nucleic acid.

8. The method of paragraph 7, wherein multiple copies of the dye bind todouble-stranded nucleic acid in direct relation to a length of thenucleic acid.

9. The method of paragraph 7 or 8, wherein the dye is an intercalatingdye.

10. The method of any of paragraphs 7 to 9, wherein the dye becomes moreluminescent when bound to the double-stranded nucleic acid.

11. The method of any of paragraphs 1 to 10, wherein the partitions aredroplets.

12. The method of any of paragraphs 1 to 11, wherein the one or morereporters include a probe that binds to the second target and not thefirst target.

13. The method of paragraph 12, wherein the data includes a signalcreated by detecting light emitted by the probe.

14. The method of any of paragraphs 1 to 13, wherein the concentrationof the first target is determined from only a first subset of the data,and wherein the concentration of the second target is determined fromonly a distinct second subset of the data or with all of the data.

15. The method of any of paragraphs 1 to 14, wherein the one or morereporters include a probe for a third target, and wherein the dataincludes data for amplification of the third target in the partitions.

16. The method of paragraph 15, wherein the step of collecting dataincludes a step of detecting, in the same optical channel, light emittedfrom the generic reporter and from the probe.

17. The method of paragraph 1, wherein all of the data represents lightemitted by the generic reporter.

18. The method of paragraph 1, wherein amplification of the first targetis detected only as light emitted by a generic reporter.

19. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing partitions containing a first target and asecond target at partial occupancy, each partition containing a genericreporter that binds to amplicons representing both targets, with aplurality of the partitions each containing a copy of both targets; (B)amplifying the targets in the partitions; (C) collecting data from thegeneric reporter for amplification of the first target and the secondtarget in the partitions; and (D) determining a concentration of thefirst target from only a subset of the data that selectively excludespartitions positive for the second target.

20. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing partitions containing a first target and asecond target at partial occupancy, each partition containing a genericreporter that binds to amplicons representing both targets and a probethat binds specifically to an amplicon representing the second target,with a plurality of the partitions each containing a copy of bothtargets; (B) amplifying the targets in the partitions; (C) collectingdata for amplification of the first target and the second target in thepartitions, at least in part by detecting light emitted by the genericreporter and the probe; and (D) determining a concentration of the firsttarget from only a subset of the data that selectively excludespartitions positive for the second target.

21. The method of paragraph 20, wherein the step of collecting dataincludes collecting data for amplification of a third target in thepartitions, further comprising a step of determining a concentration ofthe second target and the third target.

22. The method of paragraph 21, wherein the subset of the data alsoselectively excludes partitions positive for the third target.

23. The method of paragraph 21 or 22, wherein the subset of the dataincludes partitions positive for the third target.

24. The method of any of paragraphs 21 to 23, wherein each partitionincludes a probe that binds specifically to an amplicon representing thethird target.

25. The method of paragraph 24, where the probe that binds specificallyto an amplicon representing the third target does not bind to anamplicon representing the first target and does not bind to an ampliconrepresenting the second target.

26. The method of any of paragraphs 21 to 25, wherein the step ofcollecting data is performed in first and second optical channelsrepresenting different wavebands of light, wherein amplification of thethird target is detectable in the first optical channel, and whereinamplification of the second target is detectable in the second channel.

27. The method of any of paragraphs 20 to 26, wherein the probe includesa luminophore, wherein the step of collecting data is performed in firstand second optical channels representing different wavebands of emittedlight, wherein light emitted by the generic reporter is detected atleast substantially exclusively in the first optical channel, andwherein light emitted by the luminophore is detected at leastsubstantially exclusively in the second optical channel.

28. The method of any of paragraphs 20 to 27, wherein the genericreporter includes an intercalating dye.

29. The method of any of paragraphs 20 to 26 or 28, wherein the probeincludes a luminophore, and wherein the step of collecting data includesdetecting light emitted in the same waveband by the generic reporter andthe probe.

Example 7 Selected Embodiments IV

This example describes selected embodiments of multiplexed digitalassays each performed with a specific reporter and a generic reporter,presented as a series of numbered paragraphs.

1. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing partitions each including a portion of a samemixture, the mixture containing a first target and a second target andalso containing a generic reporter that is sensitive to amplification ofeither target and a specific reporter that is specifically sensitive toamplification of the second target, wherein only a first subset of thepartitions each contain at least one copy of the first target and only adistinct second subset of the partitions each contain at least one copyof the second target; (B) amplifying the first target and the secondtarget in the partitions; (C) collecting amplification data from thegeneric reporter and the specific reporter present in a plurality of thepartitions; and (D) calculating a level of each target based on theamplification data.

2. A method of performing a multiplexed digital assay, the methodcomprising: (A) providing partitions each including a portion of a samemixture, the mixture containing a first target and a second target andalso containing a generic reporter and a specific reporter, wherein onlya subset of the partitions each contain at least one copy of the firsttarget and only another subset of the partitions each contain at leastone copy of the second target; (B) amplifying the targets in thepartitions to generate a first amplicon corresponding to the firsttarget and a second amplicon corresponding to the second target, whereinthe generic reporter binds to the first amplicon and the secondamplicon, and wherein the specific reporter binds to the second ampliconand not the first amplicon; (C) collecting amplification data from aplurality of the partitions by detecting light emitted by the genericreporter and at least one luminophore of the specific reporter; and (D)calculating a level of each target based on the amplification data.

3. The method of paragraph 1 or 2, wherein each partition of a distinctthird subset of the partitions contains at least one copy of bothtargets.

4. The method of paragraph 3, wherein partitions of the third subset arenot reliably distinguishable from partitions of the first subset or thesecond subset in the amplification data.

5. The method of paragraph 3 or 4, wherein the partitions of the thirdsubset contribute to the level calculated for at least one of thetargets.

6. The method of any of paragraphs 3 to 5, wherein calculation of alevel of at least one of the targets is performed with only a portion ofthe amplification data that selectively excludes partitions of the thirdsubset.

7. The method of any of paragraphs 1 to 6, wherein all of the data iscollected from the plurality of partitions at about a same temperature.

8. The method of any of paragraphs 1 to 7, wherein the level is aconcentration.

9. The method of any of paragraphs 1 to 8, further comprising a step ofdetermining a copy number variation of one of the targets based on thelevel of each target.

10. The method of paragraph 9, wherein the partitions provided containnucleic acid from a genome, wherein one of the first and second targetsis a reference target having a known copy number in the genome and theother of the first and second targets is a target of interest having anunknown copy number in the genome.

11. The method of any of paragraphs 1 to 10, wherein the mixtureincludes primers for amplification of each target.

12. The method of any of paragraphs 1 to 11, wherein the mixturecontains a complete set of reagents to support amplification of eachtarget.

13. The method of any of paragraphs 1 to 12, wherein the specificreporter has an intact form and one or more degraded forms, wherein theone or more degraded forms are produced from the intact form duringamplification of the second target, and wherein the step of collectingamplification data includes collection of data from the one or moredegraded forms of the specific reporter.

14. The method of any of paragraphs 1 to 13, wherein the step ofcollecting amplification data is performed in at least a first opticalchannel and a second optical channel, and wherein the first opticalchannel detects light that is spectrally distinct from light detected inthe second optical channel.

15. The method of paragraph 14, wherein light emitted by the genericreporter is detected at least predominantly only in the first opticalchannel relative to the second optical channel, and wherein lightemitted by the at least one luminophore of the specific reporter isdetected at least predominantly only in the second optical channelrelative to the first optical channel.

16. The method of paragraph 14, wherein light emitted by the genericreporter is detected in each of the first optical channel and the secondoptical channel, and wherein light emitted by the at least oneluminophore of the specific reporter is detected at least predominantlyonly in the second optical channel relative to the first opticalchannel.

17. The method of any of paragraphs 1 to 16, wherein the second ampliconincludes a pair of complementary strands, and wherein the specificreporter binds to a region of only one of the complementary strands.

18. The method of any of paragraphs 1 to 17, wherein the level of theone of the targets is calculated from only a subset of the amplificationdata that selectively excludes partitions testing positive for the otherof the targets.

19. The method of any of paragraphs 1 to 18, wherein the partitions aredroplets.

20. The method of any of paragraphs 1 to 19, wherein the genericreporter includes an intercalating dye.

21. The method of any of paragraphs 1 to 20, wherein the amplificationdata is collected after an endpoint of amplification has been reachedfor each of the targets.

22. The method of any of paragraphs 1 to 21, wherein the step ofcalculating a level includes a step of plotting the amplification datato generate a scatter plot having a first axis and a second axis,wherein each axis represents a range of intensities of light detected ina different optical channel, and wherein the scatter plot includes atleast two clusters of data points with each cluster representingpartitions positive for at least one of the targets.

23. The method of paragraph 22, wherein the step of calculating a levelof each target includes a step of calculating a level of at least one ofthe targets using only a subset of the data that excludes at least oneof the clusters.

24. The method of any of paragraphs 1 to 23, wherein the first targetand the second target are not covalently linked to each other in amajority of the partitions.

25. The method of paragraph 24, wherein the first target and the secondtarget are not covalently linked to each other in the partitions.

26. The method of any of paragraphs 1 to 25, wherein all of the data iscollected from partitions at a temperature that is below a meltingtemperature of each of the amplicons.

27. The method of any of paragraphs 1 to 26, wherein the data iscollected from a plurality of partitions while copies of the firstamplicon and copies of the second amplicon in at least a subset of theplurality of partitions remain in at least partially double-strandedform.

28. The method of any of paragraphs 1 to 27, wherein all of the data iscollected with at least a portion of the first amplicon and at least aportion of the second amplicon in double-stranded form.

29. The method of any of paragraphs 1 to 28, wherein all of the data iscollected from partitions at a temperature below about 50 degreesCelsius.

30. The method of any of paragraphs 1 to 29, wherein the specificreporter includes a nucleic acid that binds specifically to the secondamplicon.

31. The method of any of paragraphs 1 to 30, further comprising a stepof classifying individual partitions as positive or classifyingindividual partitions as negative for at least one of the targets basedon the amplification data.

32. The method of paragraph 31, further comprising a step of determininga number of partitions classified as positive or classified as negativefor one of the targets, and wherein the step of calculating a level forthe one target is based on the number of partitions classified aspositive or classified as negative.

33. A composition for performing a multiplexed digital assay,comprising: (A) a plurality of droplets disposed in a same continuousphase, the droplets each including a portion of a same mixture, whereinthe mixture contains a first target, a second target, a generic reporterthat is sensitive to amplification of either target, and a specificreporter that is specifically sensitive to amplification of the secondtarget, wherein only a first subset of the plurality of droplets eachcontain at least one copy of the first target and only a distinct secondsubset of the plurality of droplets each contain at least one copy ofthe second target, and wherein the mixture includes a complete set ofreagents for amplification of each target.

34. The composition of paragraph 33, wherein the continuous phaseincludes oil.

35. The composition of paragraph 34, wherein the oil includes afluorocarbon oil, a silicone oil, or both.

36. The composition of paragraph 33, wherein the continuous phase is aliquid at room temperature.

37. The composition of paragraph 36, wherein the continuous phase is aliquid at 90 degrees Celsius.

38. The composition of paragraph 33, wherein a third subset of thedroplets each contain at least one copy of each target.

39. The composition of paragraph 33, wherein the complete set ofreagents includes dNTPs, primers for amplifying each target, and aheat-stable polymerase.

The disclosure set forth above may encompass multiple distinctinventions with independent utility. Although each of these inventionshas been disclosed in its preferred form(s), the specific embodimentsthereof as disclosed and illustrated herein are not to be considered ina limiting sense, because numerous variations are possible. The subjectmatter of the inventions includes all novel and nonobvious combinationsand subcombinations of the various elements, features, functions, and/orproperties disclosed herein. The following claims particularly point outcertain combinations and subcombinations regarded as novel andnonobvious. Inventions embodied in other combinations andsubcombinations of features, functions, elements, and/or properties maybe claimed in applications claiming priority from this or a relatedapplication. Such claims, whether directed to a different invention orto the same invention, and whether broader, narrower, equal, ordifferent in scope to the original claims, also are regarded as includedwithin the subject matter of the inventions of the present disclosure.Further, ordinal indicators, such as first, second, or third, foridentified elements are used to distinguish between the elements, and donot indicate a particular position or order of such elements, unlessotherwise specifically stated.

We claim:
 1. A method of performing a multiplexed digital assay, the method comprising: providing a mixture including a plurality of targets and reagents sufficient for amplification of each of the targets, wherein the plurality of targets includes a first target; forming partitions each including a portion of the mixture; amplifying the plurality of targets in the partitions; collecting data from a plurality of the partitions for amplification of each of the plurality of targets; identifying partition populations in the data, each partition population having a different target content of the plurality of targets; and calculating a level of the first target from only a portion of the data, wherein the portion of the data is collected from only a subset of the plurality of the partitions, wherein the subset of partitions excludes two or more of the partition populations, and wherein the two or more excluded partition populations include at least a pair of excluded partition populations that are not resolved from each other in the data.
 2. The method of claim 1, wherein the at least a pair of excluded partition populations that are not resolved from each other include a first partition population that is positive for the first target and a second partition population that is negative for the first target.
 3. The method of claim 1, wherein the step of identifying partition populations includes a step of plotting the data to form three or more clusters of data points, wherein each cluster represents at least one of the partition populations, wherein the portion of the data excludes one or more of the clusters, and wherein the one or more excluded clusters include a heterogeneous excluded cluster composed of at least two overlapping sets of data points representing partition populations having different target content from one another.
 4. The method of claim 3, wherein each data point represents one or more intensity values of light detected from each partition of the plurality of the partitions, and wherein the step of plotting the data includes a step of plotting the data with respect to at least one intensity axis.
 5. The method of claim 4, where the step of plotting the data includes a step of plotting the data with respect to at least two intensity axes.
 6. The method of claim 1, wherein the plurality of targets includes a second target, and wherein the portion of the data excludes each partition testing positive for the second target.
 7. The method of claim 6, wherein the plurality of targets includes a third target, further comprising a step of calculating a level of the second target based on only a portion of the data that excludes partitions testing positive for the third target.
 8. The method of claim 1, wherein the level is an average number of copies per partition of the first target in the subset of partitions, and wherein the average number is expected to be the same, absent statistical error, as an average number of copies per partition of the first target in the two or more excluded partition populations taken collectively.
 9. The method of claim 1, wherein the step of identifying partition populations includes a step of identifying a partition population that tests negative for each of the plurality of targets.
 10. The method of claim 1, wherein the step of identifying partition populations includes a step of identifying at least a pair of the partition populations that are not resolved from each other in a plot of the data.
 11. The method of claim 10, wherein the subset of partitions does not include the at least a pair of partition populations that are not resolved from each other in a plot of the data.
 12. The method of claim 1, wherein only a first subset of the plurality of the partitions each contain at least one copy of the first target, wherein only a distinct second subset of the plurality of the partitions each contain at least one copy of a second target of the plurality of targets, and wherein only a distinct third subset of the plurality of the partitions each contain at least one copy of the first target and the second target.
 13. The method of claim 1, further comprising a step of obtaining a first partition count and a second partition count from the subset of partitions, wherein the step of calculating a level is performed with the first partition count and the second partition count, wherein the first partition count is a total number of partitions in the subset of partitions, and wherein the second partition count is a number of partitions positive for the first target in the subset of partitions or a number of partitions negative for the first target in the subset of partitions.
 14. The method of claim 1, wherein the step of forming partitions includes a step of forming droplets.
 15. The method of claim 1, wherein the plurality of targets includes a second target, and wherein the step of collecting data includes a step of detecting light from a same reporter that is sensitive to amplification of the first target and the second target.
 16. The method of claim 15, wherein the same reporter is a generic reporter.
 17. The method of claim 1, wherein the step of collecting data includes a step of detecting light from a probe that binds specifically to only the first target of the plurality of targets.
 18. The method of claim 1, wherein the step of calculating a level is based on Poisson statistics.
 19. A method of performing a multiplexed digital assay, the method comprising: providing a mixture including a plurality of targets and reagents sufficient for amplification of each of the targets, wherein the plurality of targets includes a first target and a second target; forming partitions each including a portion of the mixture; amplifying the plurality of targets in the partitions; collecting data from a plurality of the partitions for amplification of each of the plurality of targets, the data being plottable to form three or more clusters of data points, with each data point representing a partition, and with each cluster having a different content of the plurality of targets; and calculating a level of the first target from only a portion of the data that excludes at least one of the clusters of data points, wherein the at least one excluded cluster includes a heterogeneous excluded cluster composed of overlapping sets of data points representing at least two excluded partition populations, wherein at least one of the excluded partition populations is positive for the first target, and wherein at least one of the excluded partition populations is negative for the first target.
 20. A method of performing a multiplexed digital assay, the method comprising: providing a mixture including a plurality of targets and reagents sufficient for amplification of each of the targets, wherein the plurality of targets includes a first target and a second target; forming partitions each including a portion of the mixture; amplifying the plurality of targets in the partitions; collecting data from a plurality of the partitions for amplification of the plurality of targets; and calculating a level of the first target from only a portion of the data that excludes each partition testing positive for the second target. 